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We derive a non-perturbative evolution equation for the gluon propagator in axial gauges based on 
the framework of Wetterich's formulation of the exact renormalization group. We obtain asymptotic 
solutions to this equation in the ultraviolet and infrared limits. 



*The author died in the crash of Swissair flight 111 on September 2, 1998. This manuscript was slightly edited based on a 
version circulated to a few colleagues on the day before the accident. 



I. INTRODUCTION AND SUMMARY 



Non-Abelian gauge theories, and in particular QCD, are nowadays fairly well understood in the short-distance (large- 
momentum) regime where asymptotic freedom allows reliable calcuations within perturbation theory. On the other 
end of the scale, in the long-distance (low-momentum) domain, fundamental unanswered questions remain, linked 
intimately to the phenomenon of confinement (or the lack of detailed knowledge thereof), and posing severe infrared 
problems that present a tough challenge for developing adequate non-perturbative methods to perform practical 
calculations. Whereas the non-perturbative effects on QCD Green functions are small when all relevant momenta 
are large compared to the inverse confinement length, the properties of the vacuum, the dynamics of the QCD phase 
transition, or the formation of color-neutral hadronic excitations from colored quark and gluon fluctuations, are 
completely dominated by the non-perturbative infrared physics. Although lattice simulations provide to date the 
most rigorous non-perturbative studies of QCD, they suffer in one way or another from finite lattice size effects and 
violation of translational or rotational invariance. Moreover, the continuum limit of results obtained on a discrete 
Euclidean space lattice is a difficult problem itself. 



A. Average effective action and non-perturbative evolution equation 

Therefore, it is clear that non-perturbative methods, formulated in continuous space and maintaining the symmetries 
of translations and rotations, are of fundamental need to complement insight into the infrared properties of QCD. Such 
a method has been developed |l]-[| during the last few years and has found diverse applications |J-||]. ft embodies 
the concept of the average effective action in continuous Euclidean or Minkowski space within the rcnormalization- 
group framework of quantum field theory. The basic idea is to study the theory within a volume fl oc 1 /k 4 and 
effectively integrate out all quantum fluctuations that can be localized within that volume, i.e., fluctuations with 
squared momentum q 2 larger than k 2 . The average effective action T K is formulated as a functional integral over 
the microscopic quantum fields, and can be shown to be equal to the usual effective action T for macroscopically 
averaged fields n. The vacuum properties are obtained in the limit k — ► where the volume $7 oc 1/k 4 tends to 
infinity. In this paper, however, we are interested in the non-perturbative infrared behavior of gluons propagating in 
an unconfined quarkless world. The volume of such an idealized colored world cannot, of course, be infinite, since in 
reality confinement limits it to be of the size of a hadronic state Q ~ f fm 3 . Hence, as we ignore the existence of 
the QCD phase transition between the colored and the hadronic world, we must cut out the long-distance hadronic 
physics beyond distances of order 1 fm, and need to restrict n to be larger than the mass scale of the QCD phase 
transition: 

k > K PT w 200 MeV . (1) 

As we shall see, the introduction of a new scale k into the theory is intimately related to the standard renormalization 
program of QCD, in which one needs to introduce a mass scale at which the Green functions are normalized (since 
they are not normalizable at zero momentum, due to the infrared divergence). 

The dependence of the average effective action T K on the variation of the scale k is controlled by an exact non- 
perturbative evolution equation KM, which is very sensitive to the infrared properties. It is of the generic form 



k2 8^ Tk ~ K 



K > T K 



(2) 



1 In a sense this concept is analogous to a quasi-particle picture of quantum fluctuations, wherein elementary excitations are 
effectively embodied in a quasi-particle with Compton wavelength r c oc 1/k: On distance scales r > r c the particle appears as 
an elementary object, but as one increases the resolution to shorter distances by a larger k' > k, excitations with wavelengths 
oc 1/k' reveal themselves as a substructure of the original quasi-particle. Vice versa, a decrease of resolution by lowering 
k, averages over fluctuations with longer wavelengths, and yields a larger quasi-particle. Loosely speaking, in the extreme 
short-distance limit n — > oo, the quasi-particle would be, for instance, a single elementary bare gluon, while in the opposite 
limit of inifinite volume, k — » 0, the quasi-particle would correspond to our Universe. The variation of the the scale k therefore 
controls which, and how much, physics one includes in the panorama. 



(2) 

where the kernel /C depends explicitly only on the (exact) 2-point function T K , but not on higher-order Green functions 
(which however implicitly enter in determining the 2-point function). It has been shown [0 that T K approaches the 
classical action in the ultraviolet limit k — > oo and becomes the usual effective action in the infrared limit k — > 0. 
A solution to the evolution equation (0) therefore interpolates between the (short-distance) classical action and the 

— — (n) 

(long-distance) effective action. Since T K = '^2 n T' K generates the k dependent one-particle irreducible (1PI) Green 

— (n) 

functions T K (such as the inverse propagator for n — 2, or the vertex functions for n > 3), the evolution equation 

for T K is equivalent to an infinite set of corresponding equations for the 1PI functions T K , which are the differential 
version of the well-known Dyson-Schwinger equations [7], however with an additional infrared cut-off given by k. 
Just as in the case of the infinite number of Dyson-Schwinger equations, a truncation to a finite number of coupled 
equations is unavoidable, if one wishes to find an explicit, but approximate solution. 



B. Evolution of the gluon propagator 

The purpose of this paper is to demonstrate the powerful potential of the average effective action T K and its 
evolution equation by studying the simplest non-trivial object in QCD without quarks, namely the gluon propagator. 

(2) 

Since the gluon propagator A K is related to the inverse of the 2-point function T K , we can obtain from the 
evolution equation for T K a corresponding equation for A K , which determines how the propagator changes as the scale 
k is lowered from some large initial value in the ultraviolet all the way into the deep infrared regime. Unfortunately, the 

evolution equation for A K contains in addition the unknown 3-gluon and 4-gluon vertex functions T K and T K , which 

(5) 

are themselves determined by similar, but even more complex equations, involving further higher-order functions T K , 

T K , and so forth. However, by working within the class of axial gauges, the evolution equation for the propagator 
becomes remarkably simple (at least formally), because the exact propagator is just the bare propagator times a 
renormalization function Z K , 

A K (q) = Z K (q)A^(q), (3) 

and the evolution equation (|2j) translates to an evolution equation for Z K , 

«(3) ^(4)1 



0_ 
8k 2 



K 2 ^z-\ q ) = K! ^,A K ,rr,r7 , (4) 



where the kernel K! explicitly depends on the exact propagator A K and the exact 3- and 4-gluon vertex functions 

T K and T K . In the class of axial gauges, it is furthermore possible to project out all contributions of 4-gluon 
vertex functions, so that the remaining unknown object is the exact 3-gluon vertex. The latter can be eliminated 
by exploiting the gauge symmetry properties of QCD, in particular the Slavnov- Taylor identities, which provide a 

(3) 

constraint equation between the 3-gluon vertex T K and the propagator A K . The strategy is then to construct an 

—(3) 

ansatz for T K in terms of A K such that this constraint equation is identically satisfied. As a result, one arrives at 
an evolution equation for A K in terms of the propagator alone, which must be solved as a function of k. The crucial 

(3) 

point of success in this program is the choice for T K . Although constrained by gauge symmetry, this choice is hardly 
unique. In the present paper we construct a particularly simple ansatz, since our main motivation is to illustrate the 
concept and the techniques involved. 



C. Connection of propagator with gluon distribution function 

An important point that one should bear in mind throughout is, that the gluon propagator A K (q), in general, is a 
gauge-dependent object. Only in the ultraviolet regime (q — ► oo), where asymptotic freedom is approached, it reduces 
to a gauge-independent form as given by the perturbative one-loop formula M, In the infrared domain (q — > 0), on the 
other hand, confinement should manifest itself in the behavior of the gluon propagator, but here the gauge-dependence 
foils an unambiguous assignment of confinement effects. Yet, the fact that the propagator is gauge-dependent does 
not imply that it does not contain physics; rather, it is that the physics is obscure and difficult to extract. 

Because of this problem it is desirable to relate the gluon propagator to gauge-invariant quantities, for example the 
Wilson loop or the gluon distribution function of hadrons measured in experiments. The latter is intimately connected 



with the spectral density of gluon modes in the propagator. Therefore the evolution equation for the propagator can 
be transcribed, as we shall show, into a corresponding evolution equation for the gluon distribution function. Indeed, 
in the regime where the longitudinal (or energy) component of q is much larger than the invariant q 2 , one recovers the 
famous DGLAP equation |9| , the perturbative evolution equation for the gluon distribution function. Such a physical 
scenario is realized, for example, certain hard processes occurring in high-energy hadron collisions or deeply inelastic 
lepton hadron scattering where a hard gluon can be knocked out and initiate a gluon jet with qo ~ q z ^S> q± ^S> q 2 
that evolves by means of fluctuating (real and virtual) gluonic offspring towards lower and lower momenta. 



D. Related literature 

A large body of work concerning non-perturbative analyses of the gluon propagator exists in the literature [M, 
which may be subdivided into analytical and lattice studies. 

Most analytical studies were carried out by attempting to solve the Dyson-Schwinger equation for the gluon prop- 
agator in pure SU(3) gauge theory without quarks, and in various covariant and non-covariant gauges, for example in 
the Landau gauge M-Q, the temporal and spacelike axial gauge 15 19 2l],|2(]], and the light-cone gauge J2^^j|. The 



non-covariant axial and light-cone gauges have the advantage that they are ghost-free and involve only the physical 
gluon degrees of freedom, whereas in covariant gauges one faces a complex coupling between gluon and ghost vari- 
ables. On the other hand, the structure of the propagator is more complicated in the non-covariant gauges. In either 
case, approximate solutions for the gluon propagator obtained in the literature from the Dyson-Schwinger equation 
vary widely p4] in the infrared behavior of the gluon propagator, whereas the large-momentum behavior is dictated 
by the well-known perturbative result. Predictions for the dependence of the propagator in the small- momentum 
limit include an infrared enhancement ex q~ 4 or oc q^ 2 (lnq 2 )^ 1 , infrared constant ex q 2 , or infrared vanishing ex q 4 . 
Recall however, that the gluon propagator is a gauge-dependent object, so that these very different results are not, 
necessarily, contradicting each other. 

Lattice studies are at present equally obsure, since here (in addition to the gauge-dependence) finite lattice size 
effects make it difficult to penetrate the deep infrared where the gluon wavelength becomes close to or larger than the 
linear lattice length. There have been a number of lattice simulations of the gluon propagator [p5|— p7| , all of which 
used a fixed lattice Landau gauge, and thus are plagued by Gribov ambiguities that can lead to significant systematic 
errors. It is therefore not surprising that fits to the lattice results to date are not unique and consequently do not 
allow, at present, for a definite conclusion regarding the infrared behavior of the gluon propagator. Nevertheless, 
viewed as a whole, these studies seem to suggest that the Landau-gauge gluon propagator is finite and non-zero at 
q 2 = 0, although a propagator that vanishes at q 2 — has also been claimed p5| to be consistent. 



E. Strategy of procedure 

A roadmap of our approach to arrive at a solution for the gluon propagator within the framework of the average 
effective action may be given by the following list of conceptual steps: 

1. We consider the pure SU(3) gauge theory without quarks in Minkowski space, and from the very beginning we 

choose to work in the class of axial gauges. 

2. We start from the corresponding vacuum persistence amplitude Z — exp(iW), which allows us to separate out the 

ghost contribution so that in effect we deal with a ghost-free theory involving solely the gauge fields. 

3. The generating functional W = — ihiZ is then extended to a scale-dependent version W K by including an infrared 

regulating source term 5ft K = A^VJ^ A v that is quadratic in the gauge fields A and depends on the momentum 
scale n, such that only quantum fluctuations with momenta > n are included and the limit k —> recovers the 
full theory. 

4. From W K we obtain then the corresponding scale-dependent effective action F K which generates the one-particle 

in] 

irreducible n-point functions T K , such as the inverse propagator, the 3-gluon and 4-gluon vertex functions, all 
of which explicitly depend on the cut-off scale k. 

5. Subtracting from T K the infrared regulator 3? K , and averaging over all gauge field configurations with in the 

effective volume ft oc 1/k 4 , we arrive at the average effective action F K . Differentiation of F K with respect its to 
k dependence leads the to the desired exact evolution equation. 



6. From the evolution equation for T K = ^ n T K we then project out the quadratic term T K that is related to the 

inverse gluon propagator. After decomposing the tensor structure of the inverse propagator, we obtain a set of 
coupled equations for two independent scalar functions, a K and b K . 

7. Next we focus our attention to the light-cone gauge, a special case of the axial gauges, in which the function 

a K drops out, so that we are left with a single evolution equation for the dimcnsionlcss function Z K = q 2 /b K . 
Moreover, all 4-gluon vertex contributions can be eliminated, and consequently only the 3-gluon vertex function 
survives in the determination of Z K . 

8. By constructing a specific ansatz for the 3-gluon vertex function that obeys the constraint of the Slavnov- Taylor 

identity for the gluon propagator A K , we obtain a closed equation for the A K . The formal solution of this final 
evolution equation is simply A K = Z K A K , where A^ is the bare propagator. 

9. The remaining integration of the final evolution equation for Z K must be done numerically, but in the ultraviolet 

and infrared limits, we are able to extract analytical solutions, which depend (aside from the gluon momentum 
q) on the scale k. In the limit k — > one obtains then from Z K= o(q 2 ) the full gluon propagator in the light-cone 
gauge, A(q)=Z (q 2 )A^(q). 

10. In its spectral representation, the gluon propagator can be related to the gauge-independent gluon distribution 

function G(q, k) through the renormalization function Z K (q 2 ), and the evolution equation for Z K can be tran- 
scribed into a corresponding evolution equation for G. In the high-momentum limit we recover the perturbative 
DGLAP evolution equation, and we find that our solution coincides with the perturbative result. 



F. Main results 

Although the full solution to the our evolution equation for the gluon propagator A K (q) = Z K (q 2 ) A^(q) in the 
light-cone gauge requires a numerical analysis, we are able to arrive at analytical solutions for Z K in the extreme 
limits q 2 /k 2 — ► oo and q 2 /n 2 — ► 0. In the former case, the ultraviolet limit, we obtain: 

*«(«) - x - -i&?- * UJ ■ (5) 

On the other end of the energy scale, in the infrared limit, the leading behavior turns out to be: 

v-i r 2^ « 2 ^ 9 2 Cg q 2 ,„s 

48ir z n l 
The corresponding limiting behavior of the actual gluon propagator then follows as: 

A K (q) q «°° , A K (q) 9 c7 — . (7) 

q z myq z I k ) q 

The ultraviolet behavior is consistent with asymptotic freedom, corresponding to a screening of the color charge due 
to go/g 2 = Z^ 1 < 1. The infrared solution would, on the other hand, correspond to a linearly rising potential 
V(r) ex r as r — > oo, in accordance with the phenomcnological picture of confinement. These results are certainly 
rather qualitative, firstly, because the inclusion of quark degrees of freedom which we left out here, may alter the 
details of the infrared behavior and, secondly, because the weakest point of our analysis is the aforementioned ansatz 
for the 3-gluon vertex function, which may not be all that good in the long-wavelength limit. But even for our specific 
ansatz, an exact numerical solution of the evolution equation for the propagator needs to be carried out before more 
robust conclusions can be drawn. 



G. Organization of the paper 

The reminder of the paper is structured in accordance with the above list of procedural steps: 

In Sect. II, we recall the necessary basics of the functional formalism, which we then extend to its scale(«;)-dependent 
analogue. The effective action T K for this scale-dependent functional formulation, obtained as usual, is then related 



to the average effective action T K , which is the generating functional for the Green functions in the presence of the 
cut-off k. We derive the desired exact evolution equation for the change of T K with a variation of k. 

Sect. Ill is devoted to applying the formalism to the evolution of the gluon propagator A£". We first derive, from 
the fundamental evolution equation for T K , the general equations that govern the k- variation of the propagator. Next 
we restrict ourselves to the light-cone gauge, and arrive at a considerably simpler, single evolution equation for the 
renormalization function Z K , the formal solution of which is equivalent to the solution of the gluon propagator in the 
light-cone gauge. 

In Sect. IV, we take pragmatical steps to actually solve the evolution equation, subject to a necessary assumption 
about the form of the 3-gluon vertex function. The final master equation for the renormalization function Z K and 
hence for the propagator A M „, can then be solved in closed form, and we are able to obtain the above-quoted results in 
the ultraviolet and the infrared limits. A phcnomcnological formula for the propagator that may be useful for parton 
model applications, is constructed by interpolating between the two extreme limits. 

Sect. V applies the results for renormalization function Z K to illustrate two important phenomenological connections 
with experimentally measurable quantities, namely the QCD running coupling a s (q 2 ) and the gluon distribution 
function g K {q). First, we infer from Z K the running of the coupling a s (q 2 ), using standard renormalization group 
arguments, and then we relate Z K via the spectral density p K of the gluon propagator, the gluon distribution function 
g(q, k) and its evolution equation. 

Appendix A summarizes the notation and conventions used in the paper. Appendix B recalls some basic formulae 
of the functional formalism in QCD, and provides a list of relevant Green functions and vertices. In Appendix 
C, we discuss the absence of ghosts in axial gauges, allowing a factorization of the generating functionals for the 
ghost and the gluon fields. Appendix D elaborates the details of the general structure of the gluon propagator in 
axial gauges, and the simplifications that emerge when specifically using the light-cone gauge. Appendix E briefly 
reviews the connection between the gluon propagator and the gluon spectral density, the latter being related to the 
experimentally measurable gluon distribution function. 

Table 1 provides a summary list of the notation used in this paper. 



quantity 


meaning 


A a 

T 

A » 

j; 

rf,rf 


gauge fields 

average gauge field { _4£ ) 

external gauge field current 

ghost fields 

external ghost field currents 


n 
S*[-4] 
5R«[M] 


infra-red 'cut-off 1 scale of dimension mass 

infra-red regulator that suppresses propagation of gluon modes 
with momenta q < k,. 

infra-red regulator that suppresses propagation of ghost modes 
with momenta q < k. 


Z K [J,a,a] = zi A) [J] Z^\a,a\ 
W K [J,a,a] = -tlnZ„ = W { K A) + W^ 

r„[3] = wi^U] - joA 

T K [A] = T K [A] - H„[3] 


vacuum-persistance amplitude, in presence of infra-red cut-off 
generating functional of connected Green functions Q^ 
effective action, generating functional of proper vertex functions r« 
average effective action, with infra-red regulator $t K subtracted 


/ , , \ abc 

v;l\ = (r<?>) 

/ . \abcd 
•\Aiabcd ( -pi, 4 ) 1 


exact gluon propagator (A^i'° = bare gluon propagator) 
inverse gluon propagator (alternatively: ( r« ) = I U K + II K J ) 
exact 3-gluon vertex function (V '^ v x ° = bare 3-gluon vertex) 
exact 4-gluon vertex function {W J^. c = bare 4-gluon vertex) 



TABLE I. List of basic quantities encountered in the paper. Note that all quantities with subscript k reduce to the standard 
forms when k — 0. Also, note that the separation of Z K and W K into gauge field and ghost field parts holds only in axial gauges. 



II. EFFECTIVE AVERAGE ACTION IN NON-COVARIANT GAUGES 



This Section is devoted to a brief review of the path-integral formalism for QCD in non-covariant gauges, and 
its application to the renormalization group evolution of the effective action of QCD, as developed by Reuter and 
Wetterich [g]. We refer to Appendices A and B, where our notational conventions are collected and to Table 1, which 
summarizes the notation of basic quantities encountered in the following. 

A. QCD path-integral formalism for non-covariant gauges 

We work in Minkowski space q (as opposed to the Euclidean formulation of Refs. 0), and consider pure SU(3) C 
Yang-Mills theory for N c — 3 colors in the absence of quark degrees of freedom. Our starting point is the path integral 
representation of the QCD vacuum persistence amplitude Z\J\ — (0 | Q)j in the presence of an external source J . 
Employing the conventions of Appendix B, we define the generating functional for the connected Green functions 
W[J] as usual by Z[J] = exp (i W[J}), with 



W[J] = -i\nZ[J] = -*JVln 



f VAdet(M) 5(F a [A}) exp ( * Sym [-4] + S ext [A] 



(8) 



with the normalization Af determined by the condition W[0] — pq], and 



S Y u[A] = -i ld A xT, lv T» v S ext [A] = fd 4 xj^. (9) 

Here A^ = T a A" denotes the gauge field, and the T^v = T a J 7 " the corresponding field tensor. The path-integral 
measure in (||) is short- notated as VA = n^ IIu Y\ a d^( x )t an d the gauge condition is embodied in 5 (F a [A\), 

F a [A b ll (x)} = F%(x) = 0, for all a, b, ft. (10) 

The gauge fixing determines the Jacobian det(M ) as the determinant of the Fadeev-Popov matrix 

M h (x v) - 5F ^ X "> = 5 J-A D ch 8 A (x - v) (11) 

M ab[ x,y) - ^ fc(y) SA<^ » [ Vh [ ' 

where cu b describes local gauge transformations g[u> a ] = exp (—iuj a (x) T a ), under which the gauge fields transform 
as A a — ► A^ ] a = g[uj a ] A a g^ 1 ^], so that T a T a is gauge invariant. 



Because of the practical advantages described before, we choose to work with a non-covariant gauge M,t50|, for 



which the gauge condition (10) reads, 

F%{x) - n»AXx) = 0, (12) 

where n** is a constant 4- vector, being either space-like (n 2 < 0), time-like (n 2 > 0), or light-like (n 2 = 0). The 
particular choice of the vector n^ is usually dictated by physical considerations or computational convenience, and 
distinguishes axial gauge (n 2 < 0), temporal gauge (n 2 > 0), and light-cone gauge (n 2 = 0). Among these gauges, 



2 In order to facilitate the correspondence between the functional formalism in Euclidean space of Reuter and Wetterich |2J, 
and the Minkowski space description in the present paper, we recall the translation rules between Euclidean (subscript "E") 
and Minkowski (subscript "M") formulations with metric <5 M „ = diag(— , — , — , — ) and g M „ = diag(+, — , — , — ), respectively, 

x% = (x°,x) E < — > (iar ,x) M Kb = (X,A a ) B ^^ {-iA° a ,-A a ) M 

Kb e = (D° abj r>ab) B «— » {-iD ab ,-T> ab ) M W[K,] E ^ -iW[K\ u . 

Notice that the convention for the four-potential A^, differs from that of an ordinary four- vector x^: the former is defined with 
common sign, whereas the latter has different signs of the timelike and spatial components. This is chosen for convenience in 
order to not have to change the sign of the coupling constant g when translating between Euclidean and Minkowski spaces. 



the light-cone gauge is most often employed in the literature |29). It is well suited for describing high-energy QCD 
in the infinite momentum frame | pl| , since Lorentz contraction and time-dilation causes the quantum fluctuations 
to be concentrated in close proximity of the light-cone, the direction of which naturally suggests the choice of the 
gauge vector n M . For these reasons we will later adopt the light-cone gauge by specifying n 2 = 0. For the time being, 
however, we keep n M general, so that the considerations apply to the class of non-covariant gauges as a whole. As 
elaborated in Appendix C, the gauge condition (|12|) implies for the general case of arbitrary n M 

det(M) = det(<J ac n" [3% + g ff A d J\ = det (5 ab n ■ d) , (13) 

because SF^/SA^ = S ac n M S 4 (x — y) and n ■ A^ = 0. As a consequence, the ghost degrees of freedom decouple, since 
det(M) no longer depends on the gauge field A. We may cast the generating functional (0) in a more practical form 

by rewriting the Jacobi determinant det(Af) in terms of a Gaussian integral over ghost fields rj,r], 

det(M) = jVrfDncxpli j d 4 xrj a (x) M ab n b {x)\ I VrjVn = expji j d 4 xr] a {x) (S ab n" d,,) Vb(x)\ , (14) 

and the functional 8(F%) as an exponential of a gauge- fixing action, 

VA 5(F a [A}) = VA expj-i / d A x J- (F^(x)f\ = VA expj-^ f d A x (n-A a (x)) 2 \ . (15) 

The gauge parameter £ allows here, just as in covariant gauges, to specify a particular gauge within class the of non- 
covariant gauges, e.g. Feynman-type gauges with £ = 1, or Landau-type gauges with a — 0. f\ Since det(M) in (nj) 
is independent of the gauge fields A, it can be pulled out of the functional integral over the gauge field configurations 
in QSp , so that we can factor out the ghost field dependence by rewriting (g|) as 

W[J,a,a] = -iln [z^[J] Z™[1j,ct]) . (16) 

Here, and henceforth, we have set arbitrarily the normalization Af appearing in (0) equal to unity, since it is an 
irrelevant constant factor, and introduced the functionals 

Z iA) [J] = fvAcxp(i S cS [A,J]) Z^[a,a] = f VrfDr) exp(i S eS [f),r),W,of) (17) 

with the combined gauge field action 

S cS [A,J} = S YM [A] + sfl flx [A] + S cxt [A,J] 
Syu[A}=-\ jdtxJ^FSr (18) 

sfi fix [A\ = ~ j^^^f 

S eyA [A,J\= Jd 4 xJ«AZ, (19) 

and the combined ghost field action 

S e s[rj,r],W,a] = S ghost [r), rj\ + 5 cxt [?7, ?7,ct,ct] 

«W % V]=J d A x rj a (S ab n» d„) Vb (20) 

S e ^[fj,ri,a,a] = / d*x (a a r] a + cr a rj a ) . 



3 Notice, however, that £ needs to be kept general at this point and in the following: it may be fixed only after the gluon 
propagator has been derived explicitly from inverting the terms quadratic in A in (ph. 



B. Generalization to scale-dependent formalism 

On the basis of the generating functional W[J] of (|16|), one can construct a corresponding scale- dependent functional. 
Whereas in m) quantum fluctuations with arbitrary momenta are to be included, the scale-dependent functional should 
only involve an integration over modes with momenta larger than some infrared cut-off k. A variation of k describes 
then the successive integration over fluctuations corresponding to different length scales with the aim to recover the 
full theory in the limit k — > 0. Following the rationale of Ref. S, a scale(K)-dependent generalization W K of the 
functional W in (111) is defined as 

W K [J,a,a] = W^[J] + WM[a,<r] = -i { In ( Z™ [J]) + In ( Z<"> [cf, <r]) } . (21) 



Here the scale-dependent functionals Z K are related to the usual K-independent vacuum amplitudes Z, eq. (|17|), by 
adding invariant infrared cut-offs 3? K for the gauge field A and for the ghost fields, respectively, 

Z[ A) [J] = I VA exp {z S cS [A, J] } exp { i 3? K [A] } , (22) 

Z^[a,o-]= J VffDr) expjz S eB [r),T),a,a]} exp { i 5R K [rj, n] } , (23) 

with SefilAjJ'] and Sofffry, n,a, a] defined by ( fL9| ) and (EQ), respectively, and the infrared regulators 

R«M = -lfd 4 xAZ [TZ K (d%l Al , (24) 

M K [V,V}= [d 4 xrj a [n K ((n ■ d) 2 )y b Vb . (25) 

with 



32 m° 6 cad n /q2\ / 9 M o^ \n^n v 



[K.O 2 )]™ = ^fl.0 2 )^-^ + ^J (26) 



One may wonder about the form of the infrared regulators 5ft K in ( |24| ) and (|25|). These have been constructed, so that 
they affect only the gluon and ghost propagators, respectively, as K-dependent squared mass terms that regularize 
the infrared poles in the propagators. As we shall see in detail later, by combining the quadratic pieces of SVm[»4] 
with 5ft K [-4] and similarly the quadratic terms of Sghostf??, rj\ with $t K \rj, rj\, the inverse gluon propagator A -1 and ghost 
propagator X> , respectively, are modified such that 

A" 1 ex a 2 — ► A" 1 oc d 2 R K (d 2 ) (28) 

V- 1 oc n ■ d — > V- 1 oc n ■ d (l + M^fi 'j . (29 ) 

In general, the functions R K and R K can be different, but their specific forms are unconstrained. One may therefore 
take the freedom to choose their analytic form to be the same, 

R K (d 2 ) = R K = R K , (30) 

but with different arguments d 2 , namely the operators d 2 and (n ■ d) 2 , respectively. The choice of the functional form 
for R K specifies the details of how the fluctuations with eigenvalues of the operators d 2 — d 2 and d 2 — (n ■ d) 2 larger 
than k 2 are integrated out in the computation of the path integral d2^). For example, a convenient parametrization [ 
after Fourier transformation to momentum space with d 2 — ► p -2 , and p 2 — q 2 or p 2 = (n ■ q) 2 } is (see Fig. 1), |j 



4 We shall later use a generalization of this form, which includes an additional ultraviolet cut-off A ^> k, but which contains 
© for A -> oo: R K (p 2 ) = p 2 exp (~p 2 / ' k 2 ) [exp (-p 2 /A 2 ) - exp (- P 2 /k 2 )] "\ 



R K (P 2 



P 



exp (— p 2 / k 2 



1 — exp (— p 2 /n 2 ) 
which has the following limiting behavior in the ultraviolet and the infrared, respectively: 



lim 

p 2 J K 2 — > 



R.(p 2 ) 



o 



lim R K (p 

d 2 /k 2 ^0 



(31) 



(32) 



Hence, the effect of R K (q 2 ) is vanishing in the high- momentum limit q 3> re, but provides an infrared screening as 
q — ► 0. Moreover, the original functional W of (jig), containing a// quantum fluctuations, is recovered from W K of ( |2l| ) 
in the limit re — 0, 



W„[.7,<7,<7] ^ W[y,a,o-] 



(33) 
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FIG. 1. Form of the infrared regulator i? K (g ), eq. (p2[), and its damping effect on the propagator A K (q) ex l/[g + _R K (g )] 
at small values of q 2 . The various curves illustrate the different choices of re, with k — corresponding to the case with no 
infrared cut-off at all. 



The crux of the above discussion is the convenient decoupling of the ghost degrees of freedom from the gluon degrees 
of freedom in ( pi] ) due to the choice of gauge (|J). Since we are interested in the variation of W K against re with regard 
to the physical gluon degrees, the first term in (£L]), Wk amounts to an irrelevant constant that does not affect the 

change of Wk and therefore may be absorbed in the overall normalization. In other words, for the evolution of the 
physical gluon fields with changing scale re, we can henceforth omit the ghost contribution and restrict our attention 
to Wk . Then one can derive from (plh - with reference to Appendix B - the re-dependent generalization T K of the 



standard effective action T = T k= q by introducing the average gauge field 



TJx) = (A a Jx)) K = 



sw^ 



(A) 



SJa^x) 



J=Q 



(*>(»))• 



(34) 



where the subscript re at (A) K indicates that only field modes that survive the infrared cut-off contribute to the mean 
value. The re-dependent effective action T K is then defined as the Legendre transformation of W^ 



(A) 



T K [A] = WW[J\ - Jd'xJ^A 1 : 



(35) 
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which amounts to a change of variables from { J M }, the external source, to {.A M }, the average gauge field, and yields 

5T K 



the conjugate of (|34| ) as 
8W™ 



5J£(x) 



J=o 



(aP(«o)' 



AJx) 



6a:( X ) 



A=A 



(r?» W ); 



■JZ(x) 



(36) 



Notice that switching off the external sources J = corresponds to to the extremum of the effective action at A = Aq, 
where Aq is the mean value of the gauge field for which the effective action achieves its stationary extremum [] (we 

take Aq = later). As summarized in Appendix B, repeated functional derivatives of Wk [J] with respect to the 
sources J generate the (/t-dependent) connected n-point Green functions, and functional differentiation of r K [A] with 
respect to the average fields A yields the one-particle irreducible n-point vertex functions. In particular, the second 
functional derivatives determine the 2-point functions 



-is 2 wi A) 



5j£{x)5J»{y) 



J=a 



fe 2) 

\ / f-lV 



A ab 



S 2 T K 



SA:(x)SA b (y) 



A=A 



w: - (*.->: • 



(37) 



(2) . (2) 

that is, Qk is the exact gluon propagator, and T K is its inverse, with 



K a U^v) = - 



S(Aj(y)) 

SJa^(x) 



= (A a Jx)A h M) K 



(38) 



where again the contributing field modes A are subject to the infrared cut-off at k. Similar relations hold for the 
higher n-point functions (c.f. Appendix B). 



C. Renormalization issues 

The point of introducing the scale-dependent effective action T K satisfying ( pq ) is that it allows us to vary the scale 
k, say, from some large initial value corresponding to the perturbative domain down to very small values in the non- 
perturbative regime. In effect, as we change K, more and more gluon fluctuations are included in the effective action, 
and at the same time will define the renormalizcd quantities of the effective theory, i.e., the gauge field A^ and the 
coupling g. As the effective action T K is a scalar quantity, the infinities appearing in it must take the Lorentz-invariant 
form of a scalar function Z times <S e ff[y4.] = Sym[-4] + !R K [.4], i.e., 5ym[-4] = ZSym[Ao], and !R K [.4] = Z!ft K L4o]- If we 
define the renormalized field A^ and the renormalized coupling g in terms of the bare, unrenormalized quantitities 
A 0fl and 50, 



A a 



Z 1 / 2 Al 



9o 



Z g9 



(39) 



then the bare T a ™ is renormalized by 



•M) flv — "m-^0 v O v Aq h + go J be A) LfA-Q v 



= Z 



1/2 



Ofi-A-u u v J\^ 



A/2 



WZ g )gf? c A b fi A 



(40) 



and consequently the squared field strength tensor is, 



(.F V) = z a [ {d„Al - d v A^Y + (Z^Z^gU^A^A^A: 

+ (Z A Z 2 g ) g 2 fabcfab'c'A^bA^A^A", 

Clearly, this will only take on a invariant form of a scalar function Z~ x times (J 7 ^) 2 provided that 



(41) 



5 For QCD in the absence of a medium, Aq — 0, because the vacuum is colorless and Lorentz invariant, which does not allow 
a non- vanishing (AV)- 



ff 



Z A 



7 -2 



(42) 



Indeed, as has been demonstrated originally by Kummer |32|] to order 0(g 2 ) and later been proven in general [ p3| , this 
equality of the renormalization factors for the gauge fields, the 3-gluon coupling and the 4-gluon coupling is a unique 
property of non-covariant gauges, and therefore holds in the light-cone gauge employed in this paper. Similarly, the 
infrared cut-off for the gauge fields, 3t K [-4] of eq. (g4J) is rcnormalized by 



M K [Ao] = 2< a 



Z A 



A* 



Hence, 



9^S ab R K {d 2 ) 
g^S ab R K (d 2 ) 



At 



o h 



Al 



= Z' 



sir 1 



Za 

2 



1 



-1 



(■V) 



[it-hAq b ) 


(n ■ d) 2 + R K ((n ■ d) 2 ) 1 
R K ((n-d) 2 ) 


{n v Ao b ) 




) (%X) 


\{n-df + R K ({n-df) 
R K «n-d) 2 ) 


{nuAD . 


(43) 


MM 


= Z- 1 K K [A] , 




(44) 



and consequently the forms of S c & [A] , Wh and the effective action T K [A] are preserved under simple multiplicative 
renormalization. Notice that all the physics of the renormalization group is encoded in a single scalar renormalization 
function Z, which is a function of the gluon momentum q as well as of the infrared scale k, i.e., 



Z = Z K (q) = Z K (q 2 ,n- q) , 



(45) 



where the last equality is true for the class of axial gauges, for which one can show |32] that the g-dependence can 
only enter in the combination of the two Lorentz invariants q 2 and (n • q) 2 . This function Z K will thus be the key 
to the K-evolution of the effective action and the associated gluon propagator. In particular, we shall exploit the 
advantageous property of axial gauges that (for specific choices of the gauge vector n and the gauge parameter £) the 
rcnormalized gluon propagator is simply the renormalization function Z K times the bare propagator, 



iV {q) = Z K {q) A^(q) 



1/2 

and similarly, the rcnormalized running coupling is the bare coupling constant multiplied by Z K : 

g = g(q 2 ) = Z^{q) g . 
If we choose the mass scale A as the point where we normalize the theory [c.f. ([56|)], then 

Z K (A) - 1 50 = ff(A 2 ). 



(46) 



(47) 



(48) 



The roadmap for the following is to derive a K-evolution equation for the effective action, and extract a correspon ding 
evolution equation for the renormalization function Z K , which then allows us to infer the exact propagator via (|4q ) 
the running coupling from (670, subject to the normalization condition (Eg). 



D. The average effective action 



After these preliminaries, we are now in the position to derive an average effective action T K from the effective action 
T K of (pq), as well as an exact evolution equation for this average T K within the renormalization group framework. 
This evolution equation determines how the physics changes when more and more gluon fluctuations are included in 
the functional by successively lowering k towards zero. The average effective action T K is defined [^]J|] as the effective 
action T K of (p5|) minus the infrared regulator 3? K , eq. (E4J), 



T K [A] = T K [A] - ^[A] 



(49) 



and reads in view of (p 
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T K [A\ = -i In 
= —i In 



■DA ex P { i (Syu[A] + S { £ [A] + U K [A] - ^R K [A] + S cxt [A,J] - S cxt {A,J])} 
VA cxp | i J d 4 x (-1 A£ {g^d 2 - d,d v ) A va - A (nM£) («"/") 
-^ c (^»-5^) A$A V C - \g 2 ]Tn e AlA h v A^ d 



(50) 



- ( A^V ab A ub — A h ' l? a0 ! 



W£-<) 



where 7?. K = lZ K (d ) from 129). Evaluating the functional integral, and setting A = 0, one sees that the classical 



contribution to the action S = 



5' 



YM 



S, 



(«) 



= vanishes, so that one arrives at 



A=0 



-(<>) 



■f(o) 



r K = r K [A = o] = C + r„ , 

where the T^' is the 'kinetic' part at A = 0, 

i /d 4 xA (A(°)) _1 ( a ;~ 2/ )A K (y,x), , 



p(0) 



(51) 



(52) 



and r is the 'interaction' part at A = 
1 



- g* I <Px<Fy I d 4 xid 4 yil 

2 



+ 12' 9 



'■• ■"■-WW(i l y > xi,wi)A lt (i/ 1 ,xi)A /( (i,,a:) 
d 4 xd 4 y / JJd 4 Xid 4 y, V (0) (x, Xi, x 2 ) A K (xi,yi) A K (x 2 ,y 2 ) V(y 2 ,yi,y) A K (y,x) 

+ 48 ff4 d4xd4 V Y\ d4xid4yi W^(x,xi,X2,x 3 ) A K (xi,yi) A K (x 2 ,y 2 ) A K (x 3 ,y 3 ) 

x W(j/ 3 ,j/2,2/i,y) A«(y,a;) 

d 4 xd 4 y / Y]_d 4 x i d 4: yid 4 Zi W (0) (x,x 1 ,x 2 ,x 3 ) A K (x 2 , z 2 ) A K (x 3 , z 3 ) 
J i=i 

xV(z 3 ,z 2 ,zi) A K (zi,?/i) A K (xi,y 2 ) V{y x ,y 2 ,y) A K (y,x) 



96 y 



(53) 



Here we made use of the formulae of Appendix B, in which A K denotes the exact proagator given by (|B18| ), while V and 
W are the exact proper 3-gluon and 4-gluon vertex functions given by (B20) and ( B24| ), respectively. Correspondingly, 

A« ' is the bare propagator, and V^\ W'°) the bare vertices, explicitly given by ( B22j ) and ( B26| ), respectively. The 

different contributions in T K correspond to the diagrams of Fig 2: the first term is the one-gluon loop, the second term 
is the tadpole contribution, the third term is the 2-gluon loop with exact 3-vertex, the fourth term is the three-loop 
contribution with exact 4-vertex and the last term is the three-loop contribution with two exact 3-vertices. 



Notice that the infrared regulating terms 5ft K in (|50|) affect only the contributions that are quadratic in A or A. 
Hence, if we write in analogy to (pi]) for the effective action T K in the presence of the infrared regulator 



r = r(°) j_ r 

J-k. — i-K T ■>. re , 

we have in view of (|4^) the following mapping between T K and I\: 



(54) 



(55) 
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r K = 




+ 




+ 




FIG. 2. Diagrammatic representation of the effective action F K = F K + F K , eqs. (pf|)-(j53): the first term is the 'kinetic' 
term, the second term is the tadpole contribution, the third contribution is the 2-gluon loop with exact 3-vertex, the fourth 
term is the three-loop contribution with exact 4-vertex and the last diagram is the three-loop contribution with two exact 
3-vertices. The curly lines represent the exact gluon propagator in the presence of the infrared cut-off k, the dots are 
3-gluon or 4-gluon vertices vertices, while the shaded circles (boxes) are exact 3-gluon and 4-gluon vertices. 



E. Evolution equation for the average effective action 



Following Ref. |2j], one can derive an exact evolution equation for the average effective action T K defined by (|4S|), 
which a type of renormalization group equation which governs the scale-dependence of T K as the infrared cut-off n is 
varied. Let us introduce the dimensionless evolution variable, 



' - ln( A 



cLk 1 dn 2 

dt = — = o — 
K 2 K 



(56) 



where A is some convenient mass scale at which the theory is normalized (Sees. IV and V), and which may be chosen 
to match a specific physics situation, e.g., the total invariant mass of a high-energy particle collision, or the large 
momentum transfer in a hard scattering process. Recalling (Hy), and introducing for abbrevation 



K[A,A,J] = exp i(S cB [A,J} + dl K \A}-JoA) 



(57) 



one obtains for the derivative of T K , using (|l6| - EM Ell ^||) , 



dt 



r«[A] 



dt 



wi A \j] 



- l dt ln 



VA K [A, A, J] 



[VA\A» (JfogQA, K[A,A,J] 



JVAK [A, A, J] 



J=Q 



= 2 Tr 



dt 



n, 



K ) (gP 



Ao A 



(58) 



while for the derivative of the second term on the right-hand side of (E9h one has 



^M 



JVA K[A,A,J] V 2 dt 



1 d 



ft?) ( j ' VAAy. K[A,A,J] [ VAA V K[A,A,J] 



j=o 



i T r 
2 



9 « 
dt 11 ' 



An o A L 



(59) 



where Tr [. . .] stands for the trace over all internal indices, as well as an integration over continuous variables. Sub- 
tracting ( |5S| ) from (pq), utilizing that Q k = (T K )~ 1 = (T K + TZ K )~ l , where T K is the second functional derivative 
of T K with respect to A, one arrives at the desired evolution equation for the effective average action (|4^) : 



I f ^= 2^ 



H, 



'^(2) 



. (ri ! 



K, 



= iM\ 



(60) 
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III. THE EVOLUTION EQUATION FOR THE GLUON PROPAGATOR 



Working henceforth in momentum space, we now take practical steps to solve the evolution equation (|6C|) for the 
gluon propagator. Recall that we defined the exact gluon propagator, respectively its inverse, as (c.f. (p7lT3q)): 

A K /t(g)= (G^(q,-q)) ab = (A;(q)Al(~q)) K , (61) 



(A-%l(q)= rW («,-«) 



jliJ 



r { *\ q ,-q)) +n K Z(q z ) 



(62) 



Our goal is now to infer from the general evolution equation (60) for the average effective action T K a corresponding 
evolution equation for A" 1 , from which we can then determine the properties of the propagator A K itself. 



A. The general case 



We begin by rewriting (pfl) as 

d 4 q 



! r »w - 1 



(2») 



S *"<«": 



ab 



rf(q,-q)+K K (q 2 ) 



fJLU 



ab 



= lM 



(63) 



As this is an exact equation, any attempt to solve it in full is certainly out of question, because it would require to 
solve for an infinite number of the vertex functions T K which contribute to both sides of (p3J). On the left-hand side, 
the T K enter through the series representation of T K [J5] , 



r„F] 



2-^ „\ 



d 4 q„ 



Al 



(2tt) 4 (2n) 4 



(r?W 



• Qn)) K\{qi) 



■KliQn) 



(64) 



while on the right-hand side of (|6^), the T K are implicitly encoded in the 2-point function T K . However, since we 
are here interested in the behavior of only the gluon propagator A^" = (A^A V ) K , we do not need to solve (63) for the 

average effective action r K L4] as a whole, but only for its contributions T K [A] which are second order in A on the 
left-hand side of (|63|), and which are mapped on the corresponding quadratic contributions on the right-hand side, 

.,(2 



denoted by 7^ [A] , being the second order term in the series 



iM] = £ ^ 



(65) 



That is, instead of (|63|) for the full T K , we aim at the corresponding evolution equation with respect to t 
the 2-point contributions alone, 



d_ r (2) 
dt K 



7< 2) 



In k for 



(66) 



We emphasize that (B6J) is still an exact equation: no truncations have been imposed on the way from the original 

„( 2 ) 



evolution equation (|63|). If we were to know 7k exactly, then it would be straightforward to solve for the evolution 

of T K with k. Unfortunately, the function % on the right-hand side is a tremendously complicated object, because 
it implicitly contains all sorts of contributions of higher order in the gauge fields, which one would have to determine 

by solving corresponding equations for T K , T K , and so forth. Fortunately, the gauge symmetries of QCD allow to 

relate these higher-order contributions among each other via the Slavnov- Taylor identities, and it is possible, as we 

(2) 
shall demonstrate, to obtain a closed expression for 7^ without explicit knowledge of the higher-order terms, but 

rather by their implicit inclusion through the constraint equations that follow from first principles [j. 



This is analogous to the BBGKY hierarchy [B5| of Green functions in held theory: the n-point Green functions are intimately 
coupled by an infinite set of equations of motion. For example, the 1-point function (the mean held) is determined by the 
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1. Left-hand side of the evolution equation (ffr 



Returning to (|63|), we pick out from the series representation of r K L4] in (64) the contribution T K [A] that is 
quadratic in A, and then consider A = 0, 



f^( 2 ) m 



r 



(2) 



^p-i-i/^^w^.-.)):^-.) 



(67) 



A=0 



Now, the two-point function under the integral on the right-hand side is related to the inverse gluon propagator 



"(2) 



"(2) 



by ( A k) m J = r KjM „ + H K ,^ V , since r K>flu = r^,V - ft K)/u/ and r^/nv = (3A ; ) M * = (A„)^. We may, therefore, 

(2) 

parametrize T K according to the most general tensor decomposition of the inverse propagator (A K )~J that is 
compatible with the constraining Ward identities for the class of axial gauges. This requires two independent scalar 
functions, a K (q) and b K (q), in which the g-dependence can only involve J3j] the two invariants q 2 /A 2 and n 2 q 2 / (n ■ q) 2 . 
Introducing the variable 



X 



X(n,q) = 



2 2 

n q z 



(68) 



(n ■ q) 2 
the dependence of a K and b K on q and n appears as a K (q) = a K (g 2 ,x), b K (q) = b K {q 2 ,\). Hence, we can represent 

—(2) 

r „ in the following form, 

ab 



r « («?,-?) 



/i;y 



<5 Q6 (a K (<Z 2 ,x)n-(<z) + M<Ax)Q^(<z)) , 



with the orthogonal projection operators H 

1 



Piiu(q) = 5^ + 



0^(9) = - 



i-x 



g«ft/ riMv +qu,n v n^riy 

x~H 1- x-^V- 

g z n ■ q n z 



i-x 



<7^ n^,, + g^ 



fi z 



n ■ q 



(1- x)n 2 \ n^riu 



^ 2 



(69) 

(70) 
(71) 



which have been constructed from the available vectors q^, n M and from g^ in the space n fi P ,lu = = n^Q^". In the 
absence of interactions, the bare parameters would be a K — > q 2 + R K and b K — v q 2 + R K . In general, however, the 

scalar functions a K and 6 K in (p9) embody the full information about the running of T K and, hence, of the gluon 

—(2) 
propagator which is determined by the inverse of T K , as we shall show below. 



2. Right-hand side of the evolution equation \6 



Similar as above, we need to extract from 7 K [^4] in (63) the contribution 7k [A] that is quadratic in A and then set 
A = 0. We first notice that 



Landau- Ginzburg equation, which contains the 2-point function (the propagator). The 2-point function itself is the solution of 
the Dyson-Schwinger equation, which contains the 3-point and 4-point functions. The 3-point and 4-point functions in turn 
are determined by even more complicated equations that contain higher-order Green functions. This scheme continues ad 
infinitum. The hierarchy of the equations is exact, but in order to solve it approximately, it is usually truncated to a system 
of equations involving only the 1- or 2-point functions. To achieve self-consistency of the truncated set of equations at, e.g., 
the n = 2 level, the n > 3 functions must be implicitly included by additional constraint equations. For instance, in QCD 
the Slavnov- Taylor identities relate the 3-gluon vertex function to the propagator, and can be used to eliminate the 3-point 
function. We follow such a path later in this paper. 

Here and in the following, negative powers of n ■ q axe understood in the principal value sense (gSJ, which ensures unitarity. 
Notice, that the last term oc n~ 2 in both P^„ and in Q M ^, is actually oc (n ■ q)~ 2 , as is evident from the definition of x, e q- (p8[). 



16 



7k = 7«(3 = 0] 



1 f d4q d4q ' K*( _l >\ ( d V I 1\ 



ifW; 



where Tl K (q) is the Fourier transform of (pq), 



K K a uM = R^)6 ab \ 9liV - 



i(2) 



r 



i^Gn^t^V + Q^) 



Next, we decompose T„ ' in ( |72| ) into a kinetic term (II„ ; ) and an interaction term (II re ), 

26 



(4 2) )^(9,9') = (4°) + fl K ) a (g,</) 



(72) 



(73) 



(74) 



From the relations (M), we infer 



2r(°) 



s 2 r 



S 2 ? K 



dA^SA^q') 



A=0 



A=0 



2i5 4 {q + q>) 



2i5 i (q + q r ) 



5f K . 



*A^(g) 



(75) 



Applying these to the formulae ( p0| - |53j), after Fourier transformation to momentum space, we obtain for the kinetic 



term 



ngtf («. -?) = s ab (, 2 + i? K (, 2 )) ( g, v - ^ + 1 ^ , 



(76) 



while the interaction term gives 

n '/;„,(,.-,,.= - / 7S W<£;^(«, *,-*,-*) A^(fc) 



2 J (2tt) 4 " '' 

- 2 ' d 4 fc 
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(0),acd, 



, AA' 



. dd { il\ ^)d c b [ lj 



(27TY 



V'X (l,-k,-k') Af - (k) A™ ' dd (fc') V^Ufc^-g) 



£. f d \ d \ W i0 \ acde (q,~k,-k>,-p) A^'' cc \k)Af' dd \k') 
6 J (2tt) 4 (2tt) 4 flXaT v > vj * \ j k \ / 



e a c 



d 4 k <Fp 



W ( ^ e (?,-fc,-fc / ,-p / ) A^'*'(fc) A7"-^"(fc') 



24 J (27r) 4 (27r) 4 " " A(TT 

x V S " f ',!{k,k',-p) Af>f c '(p) A^'' cd '(p') V$W>-ff) , 



(77) 



with fc' = q — fc in the second term, fc' = q — k — p in the third term, and A;' = q — k — p' , p' = q — p in the last term. 



^(2) 



(0) 



Fig. 3 depicts diagramatically the inverse propagator (T\ )^„, eq. (]7J|), in terms of the contributions II« , eq. (|52|), 
and I1 K , eq. (ph. 
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FIG. 3. Diagrammatic representation of the inverse propagator (r« ) M „ = (IT« ' + li K ) M „, eqs. Q7J)-(k7[)- The first term is 
the contribution from the 'kinetic part' IT^ , while the remaining terms arise from the 'interaction part' U K . The curly lines 
represent the exact gluon propagator in the presence of the infrared cut-off k, the dots are bare 3-gluon or 4-gluon vertices 
vertices, while the shaded circles (boxes) are exact 3-gluon and 4-gluon vertices. 



Now let us define a partial derivative <9 t * that acts only on the t = ln« dependence of LIk , but not on II K , 



a* if) 



dt ( 



nw+n, 



d? nl°) 



dt K ' 



(78) 



so that we may write the right-hand side of eq. ( pq ) in a form that is reminescent of the derivative of a one-loop 
expression, but which is exact, 



7« 



d A q 



(27T) 



- 8* 

4 °t 



in if> Z( q ,- q ) 



1 

= -Tr 
2 



? lnr l 2) 
dt 



(79) 



From this representation of 7 K we extract now the contribution which corresponds to terms quadratic in the gluon 
fields, and therefore is relevant for the evolution of the gluon propagator: We utilize 



Tr[d t *lnr( 2 )] = Tr[lnn ] + Tr [ d*( U K U^ - 1 )] - ^Tr [ 9*( U K TL™ -i U K n<?> - 1 ) ] 



(80) 



where the dots refer to higher-order terms which are cubic and higher in the gluon fields and therefore contribute 
only to the 3-point, 4-point functions, etc., but not to the gluon propagator or its inverse. Now, the first term in ( |80| ) 
amounts to an irrelevant constant which may be dropped in view of (f7q), so that we finally arrive at 



-vPO = I 

7k 2 



d 4 q 



or 



(2tt) 4 * 
1 r d A q d A p 
I J (2^)4^rj 



^l b ,(q,'q) (^ l )Z(-q,q) 



- r)* 
4 °t 



n K ,7x(q,-p) (n- l )%(-p, P ) u K %( P ,- q ) (n?)%(-q,q) 



(81) 



3. The master equations for the gluon propagator 



Now we have collected all the ingredients for the evolution equation (pq): T K appearing on the left-hand side, is 
given by d67J - f70l), and 7« on the right-hand side, is determined by (pll) together with (f76|) and (f77|). In order to 
infer from this evolution equation two independent, coupled scalar equations for the two unknown functions a K and 
b K , we project @ with P^ v and <2 M „, given by (0) and @, respectively. Using P^P* = P^, Q^xQi = Q^u, and 



P^xQi = 0, as well as n>»(I& 1 



(0) -Is 



,„ = 0, (^'(IIk^,, = 0, we obtain after some algebra, 



|a K (<Z 2 ,x) 



2 I ' //,: ' 



X "M n " 



9t 



b K (q 2 ,x) 



I — X n 
X n^n v d 



m^'-o 



X n 



dt 



nrfe-9) 



(82) 

(83) 
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We remind the reader of the complexity of these equations, which are equivalent to (f36|), and hence our comments 
after (|6q ) apply also here. The key problem becomes clear in view of (f77|), which shows that II K contains not only the 
exact propagator A K , but also the exact 3-gluon and 4-gluon vertex functions V, respectively VV. In principle, one 
would therefore have to solve even more complicated equations for V and W, and then plug the solutions into II K of 
([77]). Then (83) and (182) would contain on the right-hand sides only the unknown A K , the solution of which we are 
after. However, as we show in the next subsections, it is possible to get rid of the explicit dependence on V and W by 
(i) eliminating the 4-gluon vertices and (ii) expressing the 3-gluon vertices through the propagator A K alone. Then 
we can evaluate II K , ( p3| ) and (p2) serve to determine the functions a K and b K which, in turn, would give a unique 
solution to the exact gluon propagator from ( |74J) , 



a b 



Decomposing the propagator analogous to (|6E 



n<°> +n„ 



ab 



A 



«:£„(<?) = 5 ab (A K {q\ X ) S,M + B K (q\ X )T^(q)) 



with the projectors 



1 



X (i+^ 2 ) 



2\ Itxlv n ^v + % n v 



<t 



n ■ q 



T» v {q) = - 



i-x 



X (1 + ffl 



2 ^ %<lv n ll q l/ + q^riy 



<t 



n ■ q 



(84) 

(85) 

(86) 
(87) 



and inverting H*, + fl K on the right-hand side of (§4), one finds (c.f. Appendix D) that the functions A K and B K are 
related to a K and b K by 



A K (q 2 ,x) 



1 



a*(q 2 ,x) 



In the limit of vanishing coupling g — > 0, we have a K — > q 2 
respectively its inverse H K are, 



B K (q 2 , X ) 
R K and b K 



b K (q 2 ,x) 



R K , so that the bare propagator A K , 



A<°> (q) 



& [IV I X J- I-IV 



q 2 


+ R 

1 


M 2 ) 


Q 2 


+ R* 


(q 2 ) 



g^v 



n ■ q 



(n + iq ) — - - 



(n • <z) 2 



nl 0) ^(g,-9) == (<z 2 + i? K (? 2 )) (P^ + Q^) 



i? K (g 2 )) 



ff^ 



+ I n 2 q 2 + - I M 



e; («-<z) 2 



(89) 



(90) 



and, since A 



k = ( n« j , the following inversion property holds: A Kj m a ( TEA. + II K ) A „ = g^ v 



B. The case x 







The system of evolution equations (p3) and (p3) for the functions a K , b K) and hence for A K , B K , is still immensely 
difficult to solve, because, as is evident from (|77|), the self-energy tensor II K contains products of exact propagators A K 
(the solution of which we do not know yet) with the exact 3-gluon and 4-gluon vertex functions V and W (which are 
themselves unknown combinations of propagators). However, we can make substantial progress, if we can eliminate 
the explicit %-dependence, by considering \ = n 2 q 2 / (n ■ q) 2 — 0: There are two possibilities to achieve this condition: 
(?) choosing n 2 — > 0, or (ii) considering q 2 /(n ■ q) 2 — > 0. The first possibility corresponds to choosing, among all the 
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axial gauges with arbitrary n 2 , the light-cone gauge with n 2 = 0. The second possibility, q 2 j(n ■ q) 2 — > 0, holds for 
any n 2 , and corresponds to the quasi-real limit, by which we mean the kinematic regime in which the gluon energy qo 
is large as compared to the virtual mass \fq 2 so that the gluons are practically on-shell. Specifically, we require for 
the gluon four-momentum q^ — (qo,c[±,q z ) that 

<7o « 1 2 z » Q± » I 2 ■ (91) 

This situation is typical for high-energy particle collisions with (gluon) jet production, for example, hadronic collisions 
with center-of-mass energy E cm > 100 GeV, where the gluon (and quark) fluctuations in the colliding hadrons have 
highly boosted longitudinal momentum along the beam axis, and comparably very small transverse momentum. 
Bearing this physics picture in mind, it is then suggestive to choose the vector n^ along the preferred longitudinal 
^-direction that is dictated by the collision geometry, i.e., to choose n^ in the t-z plane, parametrized as 

Tin = (u + v, 0±,u — v) n — Auv . (92) 

The two assertions ([)]]) and (£^) imply 

q 2 

n ■ q ~ 2 v qo and — ~ . (93) 

[n ■ qY 

Consequently, from ( pq ) and (Q), we have for q 2 /(n ■ q) 2 — > or n 2 — > (assuming the functions a K and b K are finite 
for all q), 

X — 0, B ^T- — » > ( 94 ) 



so that we are left with only one unknown function A K = l/a K . We find that in this limit (]82j) and (B3J) coincide, since 

g^W = n/ = -3^^ IP", (95) 

and therefore only the tensor structure n IJ ,n v li tiU appears in both equations. Using (|95| ) together with the definition 
of Xi ec l- (|6q)i and the expression ( |77| ) for n, we obtain the master equation for o K in the limit (pjJ), 

' 2 ^ f V /tS^H^M) 



(n-q) 2 -n 2 q 2 J \ 2* J (2tt) 4 ^ 

52 j 10Y n ' {k ~ k ' )nlJ {li [ A " CT(fc) A ^>') V«^v(fc / ,fc,-g) 



2' 



where k' = q — k, and we have utilized the form of the bare 3-gluon vertex V^°\ as given by (B23|). Notice that 



(Bq) contains only the tadpole contribution and the 3-gluon vertex contribution, as diagramatically represented in 
Fig. 4: all the other 4-gluon terms that are present in Ji^ of (|77| ) vanish identically upon contraction with n^ and 
n v , because A K is orthogonal to n, which is a direct consequence of the orthogonality of II with respect to q due to 
current conservation (both properties hold, of course, also for the bare functions A^ ' and IIk ), 

n" A K , M „ = = A„, „„ n v q» 8 K , ^ = = 9 K , „ v q v . (97) 
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FIG. 4. Diagrammatics of the contraction n M n^"n„: Only the tadpole contribution (proportional to n 2 ) and the 3-gluon 
vertex contribution (proportional to ra M n„) survive. All the terms drop out upon contraction with n^ and n v , 



The initial conditions for the evolution equation (Bq) are dictated by asymptotic freedom in the ultraviolet limit as 
q — ► oo, or more precisely, q 2 —> A 2 with the normalization scale A 2 /n 2 — > oo [cf. (|5q)1: 

(98) 



which implies that the gluon propagator becomes the bare propagator at the renormalization point, 

A,(q) q2 ^ A(°)(,). 



(99) 



As we move away from the asymptotic normalization scale A, the full gluon propagator propagator ( |85| ) remains 
proportional to its bare counterpart, modulo the function A K = l/a K , which encodes all effects of including softer and 
softer gluon fluctuations in the evolution equation (|9q), 



A «, iiv{q) 



a K (q 2 ,X) 



sL(o) 



with S' given by (g6|) at the value x — 0) i- e - 



SL(l) = 3^ 



g 2 + fl K (g 2 ) 

& K (<7 2 ,X) 
n • q 



A<?> (?) , 



Hence, the bare propagator and its inverse (taking now and henceforth £ — > 0) reads 



Al°» (?) 



i? K (<7 2 



4°U«,-g) = (? 2 + iU<z 2 )) 






n^q v + q^n v ' 


n ■ q 


%<lv 




I 2 . 


) 



(100) 

(101) 

(102) 
(103) 



(0) 



and the inversion property, noted after eq. (|90|), is modified for £ — * 0: A K) p a (II* + II K ) „ = g^ — n^q v j(n ■ q). 



C. Remarks 



Let us summarize the conceptual steps of the preceding subsections. From the general form of the evolution equation 
( pq ) for the quadratic (in the average gauge field) contributions of the average effective action, we inferred a coupled 
set of equations (p5) and ( p3[ ) that determine the exact form of the gluon propagator via (J83) and ( |88| ) in terms of the 
scalar functions a K and b K . In the case of x = n 2 q 2 /(n ■ q) 2 , we could eliminate the dependence on the function b K , 
and arrived at the master equation (|9q ) for a K alone, the solution of which determines the full gluon propagator by 
simply mutiplying the bare propagator with the single function a K . The presumption x ~~ * can be achieved either by 
letting n 2 — v 0, or by considering q 2 /(n ■ q) 2 — > 0. The former possibility corresponds to going over to the light-cone 
gauge, while the latter possibility is fulfilled in the kinematic regime (VW of "quasireal" gluons. In either case, we 
have the condition (|94|), under which the master equation ( |9q) is an exact equation in the sense that it contains the 
full non-perturbative evolution associated with the function a K in general axial gauges specified by the vector n M and 
the gauge parameter £. 
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IV. SOLUTION FOR THE GLUON PROPAGATOR IN THE LIGHT-CONE GAUGE 



Recall that (pq) holds for the class of axial gauges (O) in general, viz. for any choice of n^ with finite n 2 and 
arbitrary gauge parameter £. For n 2 ^ 0, the expression on the right-hand side of this equation is then still very 
difficult to integrate, as has been discussed, e.g., in Ref. |L8| for the case PI n M = (1, 0, 0,0) n 2 = 1, and £ = 0. On the 
other hand, for n 2 = 0, which we will consider in the following, the right-hand side of (96) simplifies considerably, so 
that an exact (numerical) integration is straightforward. Moreover, we will show that it is even possible to integrate 
eq. ( Pq ) in closed form by utilizing the methods of Ref. [E2j , with the result being expressable in terms of elementary 
functions. 



A. Evolution of the renormalization function Z K for n = 



The light-cone gauge can be specified by choosing, in the parametrization (p2]), the constant vector n M , such that 
it is directed along the forward light-cone in the t — z plane. Setting in (J92j) u = and v = 1, we have 



11,, 



,,, - (1,0±,-1) n 2 = 0. (104) 

It follows then that \ = 0, and if we introduce instead of a K the dimensionless renormalization function 

„2 



z*{<f 



a,(q 2 ,0) ' 
with initial condition ( j98| ) at the normalization scale A 2 ^> k 2 in the ultraviolet: 

q 2 ^A 2 



Z.(q 2 ) 

then we may rewrite the evolution equation (Bq) as 



1 



(105) 



(106) 



1 







dt Z K (q 2 ) 



d a K (<z 2 ,0) 




dt q 2 




n»n v fd g \ 




(n-q) 2 U " ^ q 'j 




i 2 C d A k n-{k~k') v 
2 9 J (2tt)4 (n-q) 2 U 


V dt 



^(k)A K J{k')V^Ak',k,-q) 



(107) 



in which now only the 3-gluon contribution with the exact vertex function V and exact propagators A K is present, 
while the tadpole co ntrib ution, i.e., the first term on the right-hand side of (|96|), vanishes since it is proportio nal t o 
n 2 . The solution of (107) then determines the full gluon propagator in terms of Z K) so that we have instead of (100), 



A K ,^(g) = Z K (q 2 ) A(°y g ), 



(108) 



with the bare propagator A« ' given by (|102|) 



8 This case would correspond to choosing u ■ 



1/2 in eq. (B2J) 
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B. The spectral representation of propagator and vertex function 



The evolution equation (107) still contains the unknown exact 3-gluon vertex function V, which, as one would 



expect, would have to be determined first, by solving a corresponding evolution equation for V, itself involving higher- 
order vertex functions. Luckily, the gauge symmetry properties of QCD imply the Slavnov- Taylor identities, which are 
the Ward identities of QCD relating the vertex functions to the propagator. In general these relations are non-trivial, 
however, in the class of axial gauges, the Slavnov- Taylor identities have a simple form. For example, the 3-gluon 
vertex function V can be expressed in terms of the propagator A K as 



qx V ACTT (g,fc,fc') A^(fc) AY(k') = A^(fc') - K v {k) 



(109) 



where (k' = q — k). Th is Sl avnov- Taylor identity suggests the following strategy: (i) construct an ansatz for V , in 
terms of A K , such that (109) is identically satisfied, and, (ii) insert this ansatz i nto t he ev oluti on equation ( |l07| ) for 
Z K , upon which one obtains a closed equation for the propagator A K , because of ( |108| ) and (102). To do so, we adopt 
the elegant method of Delbourgo |l6| and represent the exact propagator in terms of its spectral representation 



A K ,^(q) = S' (q) 



dW 2 



P.{W 2 ) 



(q 2 + R K (q 2 ))-W^ 



(110) 



where S'(q) is defined by (101), and the singularity at W 



R K in the denominator is to be evaluated with 



the usual ie prescription. The form (110) includes the bare propagator (102), A^,^ = S 1 /{q 2 + R K (q 2 )), upon 



.(") 



setting p K (W 2 ) — 5(W 2 ). The physical interpretion of ( |110| ) is very intuitive: It expresses the propagator for a gluon 
with momentum q and subject to the infrared cut-off scale k, through the weighted spectral density p K (W 2 ) which 
corresponds to the number density of virtual gluon fluctuations with an effective mass W. The case p K (W 2 ) = 5(W 2 ) 
corresponds then to a massless, non- interacting on-shell gluon (W = 0). This notion of the spectral density p K is very 
reminescent of the gluon distribution function which is measured in lepton-hadron or hadron-hadron collisions, and 
which describes the substructure of a gluon in terms of virtual fluctuations. We will return to this issue in the next 
Section. 



Inserting the spectral representation (11C) for A K into the Slavnov- Taylor identity (109), one obtains an implicit 
equation for the 3-gluon vertex function V in terms of the spectral density p K . Since A K and p K are not known at this 
point, we must make an ansatz for V that is compatible with the Slavnov- Taylor identity. A possible form p6[ that 
satisfies the identity (109), is the following spectral ansatz: 



Af (q) A7(fc) Af (k') V X 

1 
3 



(q,k,k') 

dW 2 p K (W 2 ^ 



S^ x (q) S' v,J {k) S'P T {k') 



V^ T (q,k,k' 



[{q 2 + R K {q 2 )) - W 2 ] [(fc 2 + R K (k 2 )) - W 2 ] [(jfc' 2 + R K (k' 2 )) - W 2 



(111) 



The integrand on the left-hand side is symmetrical product of three propagators S'/[(p 2 + R K (p 2 )) — W 2 ] and the 
bare 3-gluon vertex V^ -*, weighted by the spectral density pSW 2 ). Notice that the combination of propagators and 
vertex function is just what is required to solve the identity (109), and moreover, it respects Bose symmetry, because 
all three legs are represented symmetrically. Also, the appearance of the bare vertex on the right-hand side of (111) 



does not imply that we are limiting ourselves to lowest order perturbation theory: on the contrary, the propagators 
A K attached to V^ ' are the full propagators that embody the dynami cs fr om the (perturbative) ultraviolet regime all 
the way into the (non-perturbative) infrared domain. Nevertheless, (111) is just an ansatz, and hardly unique: one 
may think of constructing a different form that is also compatible with (109) but has a richer structure R. 



9 Aitkinson et al. fll9(] have conjectured that the form (111) does not necessarily comply with the Slavnov- Taylor identity, 
because the index A of Vao-t is contracted with the q-propagator, so that it is not possible to isolate a contraction of the vertex 
function with g^. Instead a more complex ansatz is propose d in Jl9[ which avoids this asymmetry. However, in the light-cone 
gauge n 2 = 0, the ansatz of Atkinson et al. coincides with (111) for n 2 = 0, so that one may conclude that in the light-cone 
gauge these subtle ambiguities are absent. 
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C. Solution for the spectral density p K and the renormalization function Z K 



Putting the pieces together, we first multiply (107) by — (n • q) 2 = — (n ■ q) q v n u , so that both sides of the equation 
are proportional to n v '. Nex t, we multiply both sides by A K/J ^(q), in order to br ing the rig ht-h and side to the form 
A K A K A K V, as required by (111)- Finally, we insert the spectral representation (11C) and (111) for the propagators 
A K , respectivel y for A K A K A K V. As the result of these manipulations, we obtain the following equation, which 
corresponds to (107): 



d 



dt Z K (q 2 



d_ 



dW 1 



p K (W 2 



(q 2 + R K {q 2 )) - W 2 + it 



where 



V K (q 2 ,W 2 ) = q x S'^(q)U^^( q 2 ,W 2 )n" 



n ■ q 



V K (q 2 ,W 2 



Z^H' (g 2 ,^ 2 )r^, 



(112) 



(113) 



and 



n' (q 2 ,w 2 ) 



* 9 
2 



d 4 k 

WT 4 



<L aCd (9,-fc,-fc') 5' AV ' cc '(fc) S>™'- dd '(k>) V i0 Jfi C '\k',k^q) (114) 



is the self-energy function (to order g 2 ) of an intermediate virtual gluon with mass W . The remarkable feature of 
t his equation is that it is now linear in the spectral density p K of the prop agator, in contrast to the previous equation 
(107) which involved a product of propagators. After integration of (112) over dt = dn 2 /(2k 2 ) as defined by (56), the 



formal solution for Z,, is: 



1 



Z K {q 2 ) 



4°V) 



dW 2 



p K (W 2 



(q 2 + R K (q 2 ))-W 2 +ie 



V K (q 2 ,W 2 ). 



(115) 



Here the first term is determined by the initial condition (106) that Z K (q 2 ) = 1 at the normalization point A. As 
q 2 — > A 2 , the contribution Z^ ~ must reproduce the bare propagator with spectral density p K (W 2 ) — ► S(W 2 ) in the 
limit .g^O due to asymptotic freedom, i.e., 



Zi°\q 2 



dW 2 



P.(W 2 



(q 2 + R K (q 2 )) - W 2 + it 



(116) 



What remains to be done is to compute the second term in (115). Thus, we insert the explicit expressions for S' 

r {0) abc 

following expression for V K 



of (|107|), and V^J x a c of (B23) of Appendix B, into (114) for U^W 2 ^), and after some algebra, we arrive at the 



V K (q 2 ,W 2 ) = -q 2 



2ig 2 C G 



d 4 k 



n ■ (k — k') n ■ k' 



2 9 2 G G 



(2tt) 4 J ■ ' [{k 2 + R K (k 2 )) - W 2 ] [(jfc' 2 + R K {k' 2 )) - W 2 } 

= 3 re 
the integral (including a factor l/i) as I K (q 2 ,W 2 ) for later convenience. Hence, (115) becomes: 



I K (q 2 ,W 2 ), (117) 



where k! = q — k and the factor Cq = N c — 3 results from the color trace f acd f h rr1 = 8 ah Ca- We have abbreviated 



Z K (q 2 ) 



= q 



dW 2 



p K (W 2 ) 



(q 2 + R K {q 2 )) - W 2 + le 



9 2 C G 

8tt 4 



h(q 2 ,W 2 ) 



(118) 



In order to evaluate I K (q 2 ,W 2 ), we must now finally commit ourselves to a specific form of the infrared regulator 
R K (p 2 ). In general, a closed analytic solution is not possible as long as R K varies strongly with p 2 / 'n 2 , so that a 
numerical solution must be found on a computer. Specifically we would like to use a slight generalization of the form 
( |3l"| ) suggested in Sec. II, 

exp (— p z / A z ) — exp (—p^/K^) 

which includes an additional ultraviolet cut-off A»k and which contains ( pl|) for A — ► oo. Such a form introduces 

a non-linear p 2 -dependence in the denominators [{p 2 + R K (p 2 )) — W 2 ] t hat appear in (118) and (117), which 
disco urage an analytical evaluation. We inten d to investigate solutions to (118) in the near future by integrating 
(117) numerically, using the infrared regulator (|119|). 
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D. Asymptotic behavior of the gluon propagator 



Notwithstanding an exact numerical study of (118), it is desirable to obtain at least an approximate analytical 
solution in the ultraviolet and the infrared limits. This may elucidate the behavior in these two extreme limits of the 
gluon propagator A K — Z K A. K within our specific approximate approach. Furthermore, it may serve as a check for 
an exact numerical treatment. In order to extract the behavior of Z K {q 2 ) for q 2 — > and q 2 — > oo, we note that the 
dominant contribution to I^q 2 , W 2 ) of (117) arises from fluctuations a t sm all k or k' = q— k; only the presence of the 
infrared regulator R K prevents a divergence. Hence, the integrand in ( |117| ) is substantially enhanced in the infrared 
region, where k 2 ,k' 2 < k 2 , and where from (p52), R K — ► k 2 . When on the other hand k 2 ,k /2 ^> k 2 , the effect of the 
infrared regulator vanishes according to (|32|): R K — » 0. Thus, we may replace R K by 



R k (p 2 



(p = q,k,k') 



(120) 



which is independent of p 2 , as desired, but which has qualitatively the same effect as R K (p 2 ) on the propagator, in 



both the infrared and the ultraviolet, 

1 



p 2 + R K (p 2 



I/P 2 

1/k 2 



for p 2 
for p 2 



oo 




Substituting (120) in (117), we obtain 
h(q 2 ,W 2 ) « 



- / d A k 



n ■ (k — k r ) n ■ k' 



[(k 2 + k 2 )- W 2 } [(k' 2 + k 2 ))~ W 2 ] 



which can be evaluated exactly, by using the standard Feynman parametrization J3£ 
menta k = q — k' in a space of d = 2uj dimensions, 



(121) 

(122) 
and integrating over the mo- 



I^(Q 2 ,W 2 



d 2uJ k 



n-{k-k') n-k' 



[{k 2 + k 2 ) - W 2 ] \{k' 2 + k 2 )) - W 2 } 

3(n • k)(n ■ q) - 2(ra ■ k) 2 - (n ■ q) 2 



dx / cT'k 



[k 2 -2xk-q + xq 2 + n 2 - W 2 ] ' 



r(2-w) / dx(l-x)(2x-l)[ x(l-x)q 2 + k? -W 



(123) 



The remaining integral can be reduced to integrals of the type J„ dxx u 1 {\ — x) v 1 (x — y) w which are integral 
representations of the hypergeometric function F(w, u; u + V] 1/y), so that the result for I K can be cast in the following 
form I2II. 



lt\q 2 ,W 2 



(« 2 



W 2 y- 2 T{2-uj) ) 



2' ±{W 2 -n 2 ) 



1 / 3 q 2 

2 f r"' 1; 2 ; A{W 2 ~k 2 ) 



(124) 



The expression (122) is singular in d = 4 dimensions due to the pole of th e Gamma function T(2 — u>) which arises 
from the usual ultraviol et di vergence of Feynma n in tegrals of the type (122). If we wer e abl e to analytically compute 
of the original integral (117) with R K given by (119), instead of the approximate form (122) with R K replaced by k 2 , 



this divergence would be avoid ed due to the exponential suppression of momenta q > A in (119). The result (124) of 
the approximate integral ( |122] ) therefore has to be regularized by hand, which we achieve by making a subtraction at 
some mass scale /x 2 <A 2 , which we choose as /i 2 = k 2 , 



li^(q 2 ,W 2 ) = km \l^(q 2 ,W 2 ) - I^(k 2 ,k 2 ) 



(125) 



This regularized form J^ res '(g 2 , W 2 ) is then finite, because, from the following property of the imaginary part of the 
hypergeometric function F(a, /?; 7; x), 

F(a, (3; 7; x + ie) — F(a, /3; 7; x — ie) = 



2iriT(~f) 0(x-l) 
r(a)r(/3)r(l+7-a-/3) 



(x-1) 



7 — a— [3 



F(7 — a, 7 — /3; 1 + 7 — a — /3; 1 — x) 



(126) 
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one readily infers that the factor T(2 — ui) in (124) cancels in the imaginary part of th e regularized expression (125), 
while the real part is finite. Hence, the limit u> — > 2 is now well defined, and ( 125 ) can be evaluated in terms of 
elementary functions, by using some transformation properties |37J of the hypergeometric function. The result is: 



with the real part, 



l£ es \q 2 ,W 2 ) = ReI^(q 2 ,W 2 ) + ilmI^\q 2 ,W 2 ) 
Reli rcs) (q 2 ,W 2 )= - 



(127) 



(1-4*) 



3/2 



hi 



1- VI -4z 



l + Vi^I 
1 



(Az — 1) arctan . _ 
Z lV4z 



+ 3 (z ~ 1) + y ln b- 



0(1 -4s) 
0(42 - 1) 



2V3 



the imaginary part, 



where 



Im l£ cs \q 2 ,W 2 ) 



W 2 - n 2 



(1-Azf /2 0(1 -Az) 



(W 2 >k 2 ) 



Substituting (127) in equation ( |118| ) for Z K 1 , we get 

i , r „„o Pk {w 2 ) 



dW 2 



2 K (q 2 ) " y (q 2 + k 2 ) - VF 2 + ^e 1 8tt 4 



1 - 



4%/i res v,^ 2 ) 



(128) 



(129) 



(130) 



(131) 



Upon taking the discontinuity at q 2 = W 2 — k 2 , using the principal- value prescription (y ± ie) 1 = P ( - J ± mS(y), 
and calculating the imaginary part of eq. (131), one arrives at the following integral equation for p K : 



P^Q 2 



l + ^Re/^V,, 2 



Z K (q 2 ) 



9 2 C G 
8tt 3 



q 2 /i+K 2 



dW 2 



p K (W 2 



(q 2 + K 2 ) - W 2 



\mI^\q 2 ,W 2 



For the case g 2 — 0, we recover, as anticipated, the free solution for the spectral density, 



P.(q 2 ) 



% 2 ) 



2 K (q 2 ) 



= 1 



(132) 



(133) 



which corresponds to a single bare on-shell gluon. 

For the case g 2 ^ 0, we note that on the left-hand side of ( 132 ) , Relji is to be evaluated from ( 12? ) at W 2 
i.e. z — 1, while on the right-hand side of ( |132[ ) the 9- function in Im/^ from (129) cuts off the upper integration 
limit at z = 1/4, or, W 2 = q 2 /4 + k 2 . Furthermore, if we consider q 2 > n 2 (keeping in mind to let k 2 — > at the end), 



q 2 + n 2 , 



and subtract the 'single-gluon' contribution (133), to define the 'multi-gluon' contribution of virtual fluctuations, 

Ml 2 ) = P^ 2 ) - pgV), P<?V) = S(q 2 ), (134) 



we find after insertion of the expressions ( J128| — |130| ) into (132) 

H.9 2 C G 



Mi 2 ) 



i 



48tt 2 



In 



T 



9 2 C G 



48tt 2 



9 2 /4 ( 1 

dw 2 - 



'I 1 



3/2 



P^{w 2 ) , 



(135) 



where we have shifted the variable of integration on the right-hand side, W 2 



w 



W 2 ~k 2 . Notice the characteristic 



feature of the integral over w 2 : it is dominated by the contributions from the region w 2 sa q 2 , provided that p K finite 
and well-behaved in that region. From (U3q), we now can extract the asymptotic behavior of p K in the ultraviolet 



<r 



oo and the infrared q 2 — > 0. 



2G 



a) The u ltra violet limit q 2 — > A 2 (A 2 — > oo): In the large-g 2 limit, the logarithm in the brackets of the left-hand side 
of ( 135 ) dominates, so that approximately 



<7 2 Pn{q 2 . 
It is easy to see, that the form 

Mq 2 ) <= 



llg 2 C G Y j r 



48tt 2 



In (q 2 /n 2 ) 



9 2 C G 

48tt 2 



9 2 /4 



dlnw 2 w 2 p K (w 2 ) 



Ug 2 C G 

48tt 2 



(136) 



(137) 



is a consistent ultraviolet solution when substituted in (136) 



b) The infrared limit q 2 — > (k 2 — > 0): When q 2 rs k 2 with k 2 tending to kpt, eq. (ffl), we can drop the logarithm 
on the left-hand side of (135), so that, 



q 2 Ml 2 ) 



g 2 C G f q2/4 dw 2 I 1 



48tt 2 



ur 



3/2 



w 2 p K (w 2 ) 



(138) 



An approximate solution in this case is 



Mq 2 ) 



1 



CQK 



9 2 C C 

48tt 2 



(139) 



which is consistent with eq. (137) in the infrared, when 



0. 



The actual gluon propagator A K ^(q) is now obtained by inserting the spectral density (134) into the spectral repre- 
sentation (110), using the expressions for the ultraviolet limit and the infrared region, (137) and (|139|), respectively: 



Ak.^O?) 



A 



K, fllS 



(q) 



q 2 + k 2 
q 2 + k 2 



n g 2 c G (q - 

^8^ hl 



9 2 C G q 2 



for q 2 > k 2 , 



48tt 2 k 2 



for q 2 



(140) 
(141) 



where S'(q) is defined in (101). In the ultraviolet limit q 2 — ► cxd, we recover the famous logarithmic dependence 
ex 1/q 2 \n(q 2 ), while in the infrared limit q 2 — > 0, the leading behavior is a power-law ex l/q 4 . 



Th e co rresponding ultraviolet and infrared behavior of the ren orma lizatio n fu nction Z K (q 2 ) can be read off ( 14C ) 
and (141), by utilizing the relation between A Ki/i „ and Z K , eqs. (108) and ( |102| ). These asymptotic results may be 
combined into a phenomenological, but hardly unique formula which interpolates smoothly between the ultraviolet 
and the infrared limit: 



A 



q z + k z 



with 



Z K {q 2 ) = (1-C(q 2 )) 



9 2 C G q 2 
48tt 2 k 2 



C(q 2 ) 



ng 2 C G ( ,,- 

48tt 2 



(142) 



(143) 



Here the prcfactors 1 — C and C interpolate between the infrared and the ultraviolet limits, with C tending to 1 as 



q 



A 2 and approaching as q 2 — > 0, e.g. C(q 2 ) = 2[1 - A 2 /(q 2 + A 2 )]. 



In Fig. 5a, we plot this form of Z K in comparison with the asymptotic results (137) and (139), for different choices 



of k. Fig. 5b shows the corresponding gluon propagator A K in contrast to the free propagator A K 



(0) 
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FIG. 5. Left panel: The inverse renormalization function Z K , eq. (143) versus q 2 for different values of K, for different 
choices of k. Right panel: the corresponding gluon propagator A K , eq. (142), in contrast to the bare propagator A„ . 



E. Remarks 



Let us summarize the strategy that has led to the main result of this paper, name ly the asymp totic light-cone-gauge 
solutions of the ren orma lization function Z K for q 2 — > oo and q 2 — > 0, eqs. (140), (141) and ( |143[ ). We derived an 
evolution equation ( |l07| ) for Z K that involves only the exact propagator A Kj(I1/ and the exact 3-gluon vertex V^x, 
but no higher-order vertex functions. To obtain a closed equation for the gluo n propagator alone, the 3-gluon vertex 
function was related to the propagator via the Slavnov- Taylor identity ( |l09|) and constructed an ansatz for V^x, 
eq. (Ill), which obeys the constraining Slavnov- Taylor identity. The necessity of making a particular (non-unique) 
ansatz is clearly the weakest point in our approach, yet it seems to be the only way to trade in the unknown VnuX 
in order to obtain a closed equation. The resulting evolution equation (112) for Z K then contains solely the gluon 



propagator in terms of its spectral density p K , and thus expresses the intimate relation between the renormalization 
f unct ion and the full gluon propagtor (on the basis of the specific ansatz for the 3-gluon vertex). The final equation 
(112) for Z K could be solved analytically in terms of elementary functions in the asymptotic ultraviolet q 2 — > oo and 
0, provided we approximate the infrared regulator R K (q 2 ) by its asymptotic behavior in the 
-> 0, respectively. 



the deep infrared q 2 - 
limits q 2 — > oo and q 2 

The ultraviolet result (140) for Z K is characterized by the logarithmic behavior consistent with asymptotic freedom, 
since the ratio of bare and renormalized coupling constants g^/g 2 — Z^ 1 < 1, corresponding anti-screening of the 
color charge, i.e., the bare charge is larger than the renormalized one (as opposed to QED, where e^/e 2 > 1, implying 
screening of the electric charge due to virtual pair creation) . 

The infrared solution for Z K , eq. (141), on the other hand, exposes a 1/q 2 behavior, which would correspond to a 
linear static potential V(r) oc r as r — > oo, as expected for confinement in the long- wavelength limit (as opposed to 
QED, where the infrared behavior is oc 1/q 2 , corresponding to the classical Coulomb potential V(r) oc 1/r). Although 
the gluon propagator A K)(U „, and thus Z K , is a gauge-dependent object, its gauge-invariant physics content may be 
extracted by relating it to the gauge- invariant Wilson loop pq] . 



V. PHENOMENOLOGICAL APPLICATIONS 



In this Section, we apply our results to the K-dependent renormalization function Z K (q 2 ) to illustrate two important 
phcnomenological connections with experimentally measurable quantities, namely the QCD running coupling a s (q 2 ) 
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and the gluon distribution function G(q). First we infer from Z K the running of the coupling a s (q 2 ), using standard 
renormalization group arguments, and then we relate Z K via the spectral density p K of the gluon propagator, to the 
gluon distribution function G(q, n) and its evolution equation. 



A. Renormalization group equation and running coupling 

Recall that the renormalized gluon propagator, respectively the renormalized coupling satisfy, [cf. 
A K (q 2 ) = Z K (q 2 )A^(q 2 ) g(q 2 ) = Z K tf)V 2 g , 

where the scalar propagator function A K (g 2 ) is related to A K ^ u (q) by 

A^fo) = A K (q 2 ) S'^q) 



(144) 



(145) 



with S' given by (101). As in (Hq), we specify the initial conditions at the scale A>1 GeV where we normalize the 
theory, such that at q 2 — A 2 , it coincides with the bare one, 



Z K (A 2 ) = 1, 



A K ^„(A 2 ) = AW„(A 2 ), 



5o = ff(A 2 ) 



(146) 



In order to invoke the renormalization group formalism for the light-cone gauge representation of A K (q), Z K (q 2 ) 1 
and g{q 2 ), it is convenient to introduce the dimensionless propagator function D through 



^ ( JS> 5(A 2 



= ? 2 A K ( 9 2 ) 



(147) 



How the physics changes when we vary A with (k fixed) is described by the renormalization group equation for D K : 
If we change the scale A, e.g. by A 2 — ► A A 2 , then the renormalizability of the theory requires that this is equivalent 
to a rescaling of D K by the factor Z K , that is, 



/) (^, 5 (A^)j = Z K ^, 5 (A^),Aj D f^, ff (AA 2 )). 



where we have written 



zM 



3,.. [ j^, 5 (A 2 ), A 



in order to expose the implicit A-dependence in Z K . Now let us define the variable 



(148) 



(149) 



(150) 



and differentiate ( 148 ) with respect to A. Then, setting A = 1 yields the standard renormalization group equation: 



d 

OT 



D,(r,g)= 0, 



where 



Pk(9) 



d\ 



5 (AA 2 



r) K (T,g) = ^ lnZ K (^, 5 (A 2 ), A 



,2 



(151) 



(152) 



(153) 



A=l 



with j3 K denoting the Callan-Szymanzik function in the presence of the cut-off scale k, and r] K the anomalous dimension 
7] K , also being K-dependent. The solution (|151|) to the renormalization group equation for D K is obviously 



D K (r, g Q ) = D K (0, g(-r)) exp 



dr' r U (/, g(-r')) 



(154) 
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D„ 



A 2 ' 



g(A<) = D K (^,g(q<)\ exp 



q dq' 2 fq' 2 ,,* 



(155) 



which shows, since D K = q 2 A K , that the evolution of the gluon propa gator is simply governed by the multiplicative 
factor Z K involving the integrated anomalous dimension rj R . In view of (148) we therefore can make the identification, 



In Z K ( — , g{q 2 ' 



1 f dq' 2 fq' 2 



(156) 



— -l 



In order to find the large-g 2 behaviour, we return to the approximate solution Z K of (143), and invert it by 
expanding in a power series in g 2 , 



Z K (q 2 ) = 1 - 



g 2 llCg, (q 2 



(4tt) 



Inl^J + Otf 



(157) 



In the large-g 2 limit, substitution of ( |157[ ) into ( |152| ) then yields the asymptotic behavior of the /3 K -function to order 
0(g 3 ): 



PM = -/3 ( °V + 0(g 5 ) 



/?(°) 



1 nc G 

(4tt) 2 3 



The solution of (152) together with (158) then yields the (gauge-invariant) large-g 2 form of running coupling 

1 



sV) 



1 +^5 2 ln(9 2 A 2 ) 



+ 0(g*) 



1 + & I" (<?7A 2 qcd 



(158) 



(159) 



with Aq CD = k 2 exp(— l/(P^g 2 ) and g = g(n 2 )- Equivalent to (159) is the running coupling a s = g 2 /(47r) at 1-loop 
order, 



a$V) = 



12tt 



11C G In fg 2 /A 



2 
QCD 



(160) 



Similarly, the solution of (153) in the large-g limit gives the (gauge-dependent) anomalous dimension r) K (q ) to order 
0( 9 2 ): 



V(q 2 ) = 



g 2 iic G 

(4tt) 2 6 



0(g 4 



(161) 



The large-;? 2 estimates (158 - 161), resulting from our approximate solution Z K of eq. (143), agree with the standard 
results obtained within perturbation theory for the pure gauge theory p2|. 



B. Evolution of the gluon distribution function 

The gluonic substructure of a hadron can be measured in experiments, for instance in deep-inelastic lepton hadron 
scattering or high-energy hadronic collisions, through the gluon distribution function. The gluon distribution function 
is defined [|39| as the density of gluon fluctuations inside a hadron, that is, in terms of matrix elements in a hadron 
state of specific operators that count the number of gluons carrying a certain fraction x of the hadron momentum P. 
The natural choice for such a number operator would be A^A^, however, in QCD this is not a gauge-invariant object. 
Instead, one uses the gauge- invariant operator T^yT^" . The precise definition of the gluon distribution function is 
most conveniently expressed in the infinite momentum frame, in which the hadron moves in the z — t plane along the 
light cone. Employing the standard light-cone representation of four-vectors, 
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v^ = (v + ,v ,v±) 



v ± = v T = v° ± w 3 



v x = (V,u 2 ) , 



r_ = \jv1 



V V 



iV w 



v w ) — v± ■ wj_ , 



(162) 



the gluon distribution function is then the average number of gluons at light-cone time r + = in a hadron state 
\P) moving with momentum P + , with the gluon fluctuations carrying a fraction x = q + / P + in an interval dx and 
transverse momenta in a range d 2 q± [p9|: 



G(a:,q_L) 



1 



xP+ 
Here the path-ordering exponential 



dr~d 2 r± e i(q r -&**■) { P \ T +y {0., r~ , r±) £(r~,0) T\ (0,0, 0±) | P) 



(163) 



£fa,ri) = Pexpiig dr'- A+(0, r'~ , 0j_) T a 



(164) 



makes the definition (163) with the non-local operator J- +u (rj~ ) T„ {r^ ) fully gauge-invariant, as it orders the gauge- 
field operators A^T a along the line- integral between rj~ and rjj- Moreover, it provides the link to compute the gluon 
distribution function in different gauges. 

We adopt the general definition to o ur choice of light-cone gauge, for which in terms of light-cone variables the 
choice of the gauge vector n^, eq. (104) reads, 



(n+,n~,n ± ) = (0,1, Ox) 



so that the gauge constraint ( |l2|) becomes 



n ■ A = A + = A- = . 



(165) 



(166) 



Thus, the factor £ in (|163| ) is equal to unity. Futhermore, we note that specifically in the axial gauges (including the 
light-cone gauge) J^j 



T +V T+ = (d+A l ) (d+Ai) 



(167) 



where a summation over the transverse components i = 1,2 is understood, and d — d/dr . This simple relation 
involves only the transverse gauge fields Ai, which has its physics origin in the fact that in the axial gauges only the 
physical, transverse gluon degrees of freedom pro paga te, while A + vanishes and A~ is a pure gauge which decouples. 
As a consequence, the gauge-invariant definition (163) of the gluon distribution takes the following form in the light- 
cone gauge: 

G(x,qx) = xP+ I ' dr-d 2 r ± e i ( q+r ~-^- f ' ±) (P\ ^'(0,r-,f±) A(0,0,0j_) | P) , 



(168) 

summed over the transverse components i = 1, 2. In order to extend this expression to accomodate our scale-dependent 
formalism of Sec. IIB, in whic h the gluon 2-point functions carry an explicit K-dependence due to the infrared regulator 
TVj? , eq. (Pq), we generalize (168) by 






Ai U" 



n^(d 2 ) 
d 2 



A, = A* 1 



R K (d 2 

d 2 






(169) 



with R K given by (31) or (119). Thus, the K-dependent gluon distribution may be defined as 



The only non- vanishing components of the gauge-field tensor J-^" = —T vti are 



^ 



-d + A~ 



T +i = d + A\ T~ l = d-A i -d i A--ig\A~,A i 



-po = d*A J ~d J A' -ig A, A 
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G«0,q_L,<7 2 



;P+ 



l + R K (q 2 )/q 2 



dr-d 2 r ± e i(g r -&^) ( p | _4"(0, r~ , f ± ) A t (0, 0, 0_l) | P) 



(170) 



Now, as discussed in Appendix E, the expectation value of the gluon number operat or A l Aj on the right-hand side 
is essentially the gluon spectral density p K that enters the spectral representation (11C) of the gluon propagator. 



Precisely, it is the transverse spatial component p£~i = Pk + Pk , of the causal correlation function 

Pk + \M = /<**re* r (P\A»(r)A»(0)\P). 
at r + = 0. Similarly, the anti-causal correlator is defined as 

P <jrUq) = \d\e^ (P\A"(0)A»(r)\P). 



The spectral density is the sum of both contributions, 

1 



Pn,^(<l) = 2 [P^\M) + Pl,l»(l)\ = P^\u{l) 



(171) 



(172) 



(173) 



where the latter equality holds only if translational invariance is preserved (in which case the crossing relation p K , ^„ = 
— pk, Ifi — pk, \iv exists), while it is invalid in physics situations where one encounters a spatially inhomogenous medium. 
In the prese nt co ntext, we are interested in the gluon distribution of a physical hadronic state in free space, so that 
we may use ( |173| ) to relate the spectral density to the gluon distribution. To do so, we first note that in the light-cone 
gauge, the tensor stucture of p K)Ml , is identical to that of the propagator [cf. Appendix E], 



p K ,»v(q) = p K (q 2 ) 9tLv 



Defining the density p K (x, qj_, q 2 ) through 



n^qu + q^n v 
n ■ q 



Pk(q 2 ) S'{q) 



P.(q 2 ) 



dxd 2 q± p K (x,qj_,q 2 ) , 



(174) 



(175) 



we see from (17C - 175) that the spectral density p K (x,q±,q 2 ) can be identified with the gluon distribution (170), 



p K (x,q±,q 2 ) = .g K (a;,q_L,g 2 ) , 



(176) 



as one may intuitively expect, since the gluon distribution measures the density of gluonic fluctuations which is nothing 
else but the 'level density' of gluon states described by the spectral density. Accordingly, the bare density pk , in the 
absence of interactions, just 



-,(0) 



p^{x,^,q 2 ) = gW(x,q^q 2 ) = 6(1 - x)5 2 (q ± ) S(q 2 ) 



(177) 



corresponds to single bare gluon carrying the full momentum fraction x = 1. We remark that the density (17£) 
satisfies the following sum rule [gjj , 



/ dxxp K (x,q ±7 q 2 ) = 1 
Jo 



(178) 



As an immediate consequence of the above identification of p K with the gluon distribution g K , the evolution of the 
latter is governed again by the renormalization function Z K : Since the gluon propagator A K ^„(q) = Z K (q 2 )Ak, )J ,v(q), 
we see from (110) that also p K (q 2 ) — Z K (q 2 )pk (q 2 ), To derive the precise form of the evolution equation for g K , let 



us consider the transverse-momentum integrated density, 



Pn{x,q 2 ) 



\ d 2 q± p K (x,q±,q 2 



(179) 



and introduce the ^-moments 
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MN,q 2 ) 



dx x p K (x,q ) 



The first moment is just 

p K {q 2 ) = p K (hq 2 ) = Z K (q 2 )^(l,q 2 ) 
Expressing Z K in terms of the anomalous dimension rj K , eq. (|156|), 



Z.(q 2 ) 



\f d 4^n«{q'\9{q> 2 )) 

1 Jk 2 q 



the iV-th moment generalization of (181) may be written as 
p K {N,q 2 ) = exp 



\[ 2 ^r^( N ^' 2 '9(q' 2 ))\ P^HN,q 2 ) = Z K {N,q 2 )^{N,q 2 ). 
The evolution with q 2 of the spectral density in TV-space is therefore governed by the evolution equation 



<) 



1 



dq 2 r^"'* ' 2 

If we define a probability distribution P (x, g(q 2 )) as 

dxx N - l P{N : g{q 2 )) 



q 2 —p R (N,q 2 ) = -~ VK (N,q 2 ,g{q 2 )) p K {N,q 2 ) 



■7 lK (N,q 2 ,g(q 2 )) , 



we can express (184) as follows 



r L d z 

q 2 g-^ Pn(x,q 2 ) = / — P(x,g(q 2 )) p K (x,q 2 ) 



(180) 



(181) 



(182) 



(183) 



(184) 



(185) 



(186) 



This evolution equation has the form of the DGLAP master equation M, however, with the essential difference that 
it contains the non-perturbative infrared physics as well, while the DGLAP equation corresponds to the perturbative 
limit of (186). This is easily realized by expanding the probability function P in power of g 2 , 



P(*,g(q 2 )) 



g 2 (g 2 ) 

8tt 2 



p{o){x) + i^W) 2 p{l){x) + 



(187) 



and substituting in ( |186| ). It is now evident that Pw must coincide with the DGLAP probability for gluon splitting, 
9~> 99 §], 



P(°)(a 



P 



DGLAP 



g^gg 



(x) = 2 C G 



1-x 



1 — X- 



+ z(l-x) 



(188) 



Hence, one may regard pQ| P (x,g{q 2 )) as a generalization of the DGLAP probability to all orders in g 2 /8ir 2 , or 
a s /2ir. 



The integral form of the evolution equation (186) can now be expressed as 

p K (x,q 2 ) = Z K (q 2 ) P { ° ] M + Z K (q 2 ) £ ^ £ | P ( z , g (q> 2 )) p K (£,g») Z? (^) , (189) 



where p K is defined in (177). Multiplying by x and integrating over x from to 1 yields on account of the sum rule 
(|178|) an integral equation for Z K in terms of the probability P, 



Z~\q 2 ) = 1 



dq' 



2 q' 2 jo 



dzP{z,g{q' 2 )) Z-^{q 2 ) 



(190) 



This equation is reminescent of (118) encountered in the context of the evolution of the gluon propagator, reflecting 
the universal role of the renormalization function Z K in the light-cone gauge. 



33 



ACKNOWLEDGMENTS 

This work was supported in part by the U.S. Department of Energy under contract DE-AC02-98CH10886. 

APPENDIX A: DEFINITIONS AND NOTATION 

This Appendix gives a summary of the basic quantities encountered in the paper, and the various notations used. 
Throughout the paper pure SU(3) C Yang-Mills theory in Minkowski space is considered, with N c — 3 colors and the 
absence of quark degrees of freedom. 

Our convention for placing indices and labels are the following: 

• Lorentz vector indices /i, v, . . . may be raised or lowered according to the Minkowski metric g^ = 
diag(l, — 1, — 1, —1), and the usual convention for summation over repeated indices is understood. 

• Similarly, color ind ices a, b, . . . may be raised or lowered according to the commutation rules of the SU(3) 
generators, eq. (jAq). 

• All other labels that do not refer to internal degrees of freedom, as e.g., T K or T^ 2 \ are consistently placed either 
as subscripts or superscripts. 

In order to avoid 'inflationary labeling' with sub- or superscripts, we often choose to suppress the color indices of 
vectors or tensors, when the color labels corresponds to the Lorentz indices, e.g., T°^ v (q, q') = T^. 

Furthermore, the following shorthand notations are employed: 

A-B = A IM g>"'B 1/ , AoB = fd i xA l ,{x) B fM {x) (Al) 

(AB)-C= A»B V C„ V , (AB)oC = J d^xd^y (A»(x)B v (y)} C^(x,y) . (A2) 

We use the symbol Tr for the trace over discrete indices, 

( A a ^(x)B b l ,(y)gf n '6 ab S' i (x - y) in space-time 

Tr [AB] = i . (A3) 

[ Af l (k)B b (k')g^S ab 5 4 (k - k 1 ) in momentum space 

Similarly, we use the symbol Sp for tracing over discrete indices as well as integrating over continuous variables, 

{/ d 4 xd 4 y Tr \A(x)B{y)] in space-time 

. , , • ( A4 ) 

/ IJ^AfF Tr [A{k)B(k')] in momentum space 

The gauge field is denoted by A^ (x) = T a A't (x) , and the corresponding gauge field tensor and covariant derivative 
are defined as 



1 

Hs) 

D„(s) = B„ - igT a Al{x) =d„- igA^x) , (A5) 



F lu ,{x) = T a J*,(x) = -—-[D,,D 



or, explicitly in color components, 

D a ; = 8 ab d,-gr c b A^. (A6) 



ib rata „ rab 
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The derivative d^ = d/dx^ acts on the space-time argument x^ — (x°,x), and the generators of the SU(3) color 
group are the traceless hermitian matrices T a with the structure constants / , as matrix elements (a, b, . . . running 
from 1 to N c ) with 

Tr(T a ,T 6 ) = N c S ab , [T a ,T b ] = +i f abc T c , (T a ) bc = -»/£ = -i f abc . (A7) 

For example, the Yang-Mills action reads then with these conventions: 

S YM = -\Jd 4 xF« u (x)F^ a (x) 

d 4 x[ {d^Alf - {d li A a v ){d v A u ' a ) + gfabc{d^Al)A^ b A^ c + g 2 f abc f ab ' c ' A^AIA^' A^ c ') (A8) 



APPENDIX B: SCALE-DEPENDENT GENERATING FUNCTIONALS AND TV-POINT FUNCTIONS 



Here we recollect the formulae for the various functionals, Green functions and vertex functions that we refer to in 
the paper. We restrict ourselves to the case of non-covariant gauges and focus our attention on the gauge field sector. 
Our formulation is in complete analogy with the usual pathintegral formalism of QCD, except for the presence of the 
infrared scale K which effectively truncates the theory to one which includes only field modes with momenta > k. 
In the limit k ^ the full quantum theory is recovered, whereas the opposite limit k — ► oo correponds to the pure 
classical Yang-Mills theory. 

The scale-dependent vacuum persistance amplitude Z K \J\ = (0 | )j t K in the presence of an external source J and 
the infrared regulator 3? K (with lim K ^o 5Rk = 0) is defined as, 



Z K [J] = N' / XUdet(M) 5{F[A}) exp $ [S Y m[A] + Jo A) 



exp 



(*K«H]) , 



(Bl) 



and the expectation values of time-ordered products of field operators (in the presence of 5ft K ) are given by, 
(A%(x 1 )...A'Zr n (x n )) K = (0\T[A;\(x 1 )...A«:(x n )]\0) ti 



(B2) 



-^ - j VA det (M) S (F[A]) exp [ i (S YM [A] + J o A)] exp ( * X K [A]) T [A% (si) . . . Afc (x n ) ] 



Z K 



Here the functional integration is over all gauge field configurations with the path-integral measure VA = 
Y\ x Y\nY\ a dA^( x ): an d Sym[-A] = —\\d ii xT il , v T lll/ . The determinant det(M) is the Fadeev-Popov determinant for 
the matrix M a b(x,y) — SF^(x)/Sco b (y) with the gauge constraint for non-covariant gauges F a [A{x)\ — n ■ A a (x) = 
(n^ being a constant 4- vector). As discussed in Sect. 2, the factor det(M) 5 (F[A)) can be converted into a ghost field 
contribution to the action in the exponential of (Bl). The grea t ad vantage of non-covariant gauges is the decoupling 
of the ghost degrees of freedom from the gauge field, so that (Bl) can be written as a sum of a ghost contribution 
and a gauge field contribution, 

Z K [J,a,cr] = Z^[J] + Z^[a,a], 



where [ c.f. @-@> ], 

Zi A) [J] = / VA exp 
zW[a,a] = f VA exp 



d"x 



-770 ■jrfj'V 



^ (n»A«) 2 



d 4 x(rj a (5 ab n»d^) n b + a a r, a + o- a rf 



j;;Ai:) \ exp(iR„[-4 



exp 



(B3) 

(B4) 
(B5) 



where 5ft K ^4] and 5R K [?7, rj\ are given by (pj) and (pa), respectively. Concerning the dynamics of the gluon gauge fields, 
the ghost contribution amounts to a constant term that factors out when generating the gluon Green functions from 
(B3) via repeated functional differentiation Z^^O] 8 n Z K [J,a,a)/8J n \ T_^_„_ n - F° r the same reason, the normal- 



J=a=a=0' 



7{A) 



ization Af' in (Bl) is irrelevant. Hence we focus on the pure gauge field functional Z\ , eq. (B4), and define for 
convenience 



S eS [A,J] = Jd*x 



1 I 9 

f a fH v — (r,V A a \ 

4 ^txv^a 2£ V A M/ 



qa Ail 



(B6) 
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1. The functional Z K [J] 

We write the gauge-field part of the scale-dependent vacuum persistence amplitude as, 

Z {A) [J] = fvAcxpU (S e &[A,J]+$t K [A]) 



(B7) 



The gluon n-point Green functions Q K (including both connected and disconnected parts) are then defined as the 
expectation values of time-ordered (T{. . .}) products of n gauge fields in the presence of the infrared regulator 3? K , 



(0l n) ( 



■n <;,)) = ( A a ^( Xl )...A a ^(x n ) ) 

fll—Hn 



s n zi A) [j] 



HT 

Z (A) [J] <*X" MW£-i (*n-l) ■ ■ ■ 8J£ (an) 

1 



J=0 



z (A) [j] 



VA exp 



i (S cB [A,J]+$t K [A]) T{A%{x 1 )...All{x n )} 



J=o 



such that the Volterra series representation of Z reads 



»=o 



Z (A) [J] = J2 h d ^n-..d A Xl (gjr>(x 1 ,...x n j) 1 '""* ,J a ^(x,)...J a ^(x n ). 

— ' m j \ //ji.../j„ 



(B8) 



(B9) 



2. The functional W K [J] 



Corresponding to (B7), we define the scale-dependent connected Green functional as, 



Wi A) [J] 



\nZ^[J] 



i In 



VA exp 



i (SMtf+ft^A]) 



W K generates connected n-point Green functions Q K in the presence of the infrared regulator !R re , 

a± ...a n 



(gv>(xi,...xn) 



( A%(xi)...A%(x n ) ) 



(c) 



= H)"" 1 



5 n W {A) [J] 



J=0 



SJ£ {x n )8j£? (*n-l) ■ ■ • 8J£ (si) 

<$X" (l.)*-! (*n-l) • • • W M V ^ [ Jl SJ £ ^ ) , 



(-?) 5zr j [j] 



(-A)r 



.7=0 



which generate the Volterra series 



W (A) [J] = ]T % —r- d i X n ...d i X 1 (gW(xi,...X n )) ai "' ,J^l)...J^(Zn). 



(BIO) 



(Bll) 



(B12) 



3. The effective action r K [A] and average effective action r K [.A] 

The effective action is the generating functional for the proper vertex functions. It is obtained as usual from Legendre 
transformation of (BIO), by defining the average gauge field A (as opposed the the gauge field A of the quantum 

fluctuations), A (x) = SWk [J]/8J£(x) = (A a ^(x)). Then, the transformation of W K yields the scale-dependent 
vertex functional T K , i.e., the effective action in the presence of the infrared regulator 5ft K , 
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T K [A] = W^ - Jo A 



i In 



VA exp 



(S e s[A,J}+?H K [A)) - Jo A 



(B13) 



One may think of (B13) as a change of variables from {J(x)} to |v4(x)}, the latter being the natural variable 
of the Legendre tranform r«[A|. The derivative of r re [A| with respect to A gives the Legendre conjugate relation 
5T K [A]/6A a (x) = —Ju(x). Repeated functional derivatives of r K [A] generate the one-particle irreducible rt-point 
functions, or proper vertices, at the stationary point A = Aq that maximizes the effective action 1\L4], corresponding 
to vanishing sources J = 0: 



S n T K [A] 



SA al {Xn)5Aal_l (x n -l) ■ ■ ■ SA^ (x{) 



(rl n) ( 



X\ , . . . X n ) 



A=A 



a\ ...a n 



(B14) 



The series representation for T K reads then: 

r M = E -1 d*x n ...d A Xl (ri n \x 1 ,...x n )Y 1 l X 1 M)---Xl^n) 



(B15) 



n=0 



Finally, the average effective action T K is defined as the effective action T K of ( B13 ) minus the infrared regulator 5ft K 
at A, 



T K [A] = T K [A] - M K [A] =-iln VA exp i (S cS [A,J] - J o A + to K [A]) - 5R K [A 



(B16) 



4. n-point Green functions and proper vertices for n < 4 



Using the above definitions of the generating functional Z K , W K , T K and T K , we list below the associated Green 
nctions G^\ g<?\ Tt\ a. 

The 1-point functions read: 



functions gi n \ Q& , T^\ and rj° for n = 1 .. .4. 






(B17) 



The 2-point functions are given by: 



r a , , — r b , 



<^(z,y) =( A;(x)Al(y) ) K = Af^y) + ^(x)^(y) 



(rl 2 )(x,j/))°' =[A- 1 ]^(x, y ) 



where the exact gluon propagator (% and its inverse T K , are defined, respectively, as 



SJa^(x) 



S{AZ(x)) K 



(B18) 



(B19) 
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For the 3-point functions one obtains: 

(& 3) (*,y,*))°^ = ( A«(x)Al(y)A c x (z) ) K = (_*)^|^ Aftfc,,*) + ( A^(x,y) (^(z)> K + 



g?Hx,y,z)) = ( A^x)Al(y)A c x (z) )^ = (-»)_£_ A^fo,*) 

> V.v\ F 0ja(X) 

(rW(x,y,z)) abC = V;l\{x,y,z) , (B20) 

where the function V is the exact proper 3-gluon vertex, 

-igV$,(x,y,z) = -igV f; v (x,y,z) + 0(g 3 ) , (B21) 

which, to lowest order in the coupling constant <?, reduces to the bare 3-gluon vertex Vq, 

V f^(x,y,z) = f abc {gx^d y -d x ) v 5\x,z)5\y,z) + g„ v (d z - d y ) x 5 i (y,x)5 4 (z,x) 

+ g„x(d x -d z )^5 4 (x,y)S\z,y)} . (B22) 

In momentum space, it reads, 

V l%{k u k 2 , k 3 ) = -i f abc { g^ih - h)„ + g^(k 2 - k 3 )x + g„x(k 3 - k 2 )^ (B23) 

Finally, the 4-point functions have the following forms: 

(g^(x,y,z,w)) abCd =(A;(x)A b MA x (z)A d a ( W )) K = {-if *^- , A£(z,w) + 

+ H) ^R ( A "V( S ^)^W) K + A b ^ x (y,z)(Al(x)) K + 
+ A c K X(z,x){Al(y)) K ) + {Al(x)) K {Al{y))M^))M d A^)). 
(gW(x,y,z,w)Y bCd = < A^x)A b v (y)A c x (z)Ai(w) )^ = {-i)*—*L—A£{z,w) 

(r^(x,y,z,w)) abCd = W$&(x,y,z,w) , (B24) 

with the function W denoting the exact proper 4-gluon vertex, 

-g 2 Wl% d Ax,y,z,w) = - g 2 W f^(x,y,z,w) + 0(g 4 ) , (B25) 

which, to lowest order in the coupling, defines the usual bare 4-gluon vertex Wo, 

w frf a (x, y ,z,w) = - { (f ace f bde -f ade f cbe ) gx^« + (f abe .f cde -f ade f bce ) gx.g^ 

+ (r e f dbe - /«*/«*) gx<ygv ^ S 4 (x, y)S 4 (z, w)6 4 (y, z) . (B26) 

In momentum space, it reads, 

W f^(ki,k2,h,k 4 ) = - { (r ce f bde -f ade f cbe )gx^a + (f abe f cde -f ade f bce )gx,g^ 

+ (f ace f dbe -f abe f cde ) gx*g^} ■ (B27) 
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APPENDIX C: FADEEV-POPOV DETERMINANT AND DECOUPLING OF GHOSTS IN THE 

LIGHT-CONE GAUGE 



In this Appendix the standard procedure of gauge field quantization is applied to the class of non-covariant gauges 
(O), and it is shown that ghost degrees of freedom are indeed absent, reducing the general non- linear dynamics 
in of QCD essentially to a linear QED type dynamics. We mention that an alternative, non-standard method was 
originally suggested and carried out in detail by Kummer B2], which elegantly avoids the ghosts altogether and 
instead introduces a Lagrange multiplier field that carries the fictious degrees of freedom. For an excellent review and 
bibliography, see Ref. E9] . Recall that under local gauge transformations 



the gauge fields transform as 



g[0 a ] = exp(ig0 a (x)T a ) , 



A a _^ A (0)a = g^ A a g -l [0 a h 



(CI) 



(C2) 



implying that Pp V J-£ v — T^J a 3-\iJ a , and thereby ensuring the gauge invariance of the Yang-Mills action Sym [-4] = 
— j J d^xT^vT^. However, a source term of the form J o A is not gauge invariant under the transformations (|Cl|). 
Consequently, the functional 



zi" 



= / VA exp i (S Y m[A] +JoA) 



exp 



(*»44) 



(C3) 



is also not a gauge invariant quantity. As is well known, this can be remedied by applying the formal Fadeev-Popov 
EJ procedure and integrate in the path-integral Z K over all possible gauge transformations g(9 a ) subject to the linear 
subsidiary condition 



F[AW] = n»A^ a (x) - (3 a (x) 







(C4) 



with normalized space-like vector n M and (3 a (x) an arbitrary weight function. The Fadeev-Popov trick to implement 
the constraint (C4) in the non- invariant functional Zk 1VB by multiplying with 

Ve]Js[r[A^]) det(M), (C5) 

a 

where the determinant is the Jacobian for the change of variables 4> a — * Q a 1 

b a [A^ } }\ 



det (M ab ) = dct 



^[A^}=0 



ve\[5U a [A^] 



Following this procedure one arrives at 

Z K = f VAdet(M) J[6 (fj> a [Ap]) exp[ i (S YM [A} + J ° A)] exp ( iU K [A] 



(C6) 



(C7) 



which is now a gauge invariant expression due to the proper account of the subsidiary condition (C4) that guarantees 
the correct transformation properties of the gauge fields in the presence of the sources J . 

To obtain the final form o f Z K as quoted in (16), one carries out the functional integration over the arbitrary 
fruitions j3 a (x) introduced in (C4), by choosing, e.g., a Gaussian weight functional 



w[d a ] 



exp 



~ / d'x [/?<V l2 



(C8) 



with the real-valued parameter £, upon which the Fadeev-Popov determinant det(M) can be rewritten in a more 
suitable way: 



det(M) = JV(3 TJ expj-^ J d 4 x [f3 a (x)] 2 ^ 6{n^A^ a (x) - p a {x) 



(C9) 
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In order to calculate the determinant, it is sufficient to integrate over 9 a in a small vicinity where the argument of 
the 5-function passes through zero at given A^ 6 ' a and (3 a . For infinitesimal gauge transformations 

g[9 a ] — > Sg[6 a ] = 1 + ig6 a {x)T a , (CIO) 

the gauge fields transform as 

K —> K + SA% , 5AI = ff/fc^ + d»e a , (Cll) 

so that one obtains 

i(n">lf 8 (i) - /? a (x)) = {(nMf°(i) + gf% c d b n»A { ° )c + n^d^6 a - /3 a ) 

= S(gf^9 b (3 c + »"^fl a ) , (C12) 

because n^Ajj. — (3 a . This latter expression is evidently independent of A't. Therefore, when substituted into (CS) 



and the integrations carried out, 

det(M) = det ( <S ac ra" [S b a d„ + g f? A*]\ = det (S ab n ■ d) , (C13) 

one sees that det(M) is also independent of the gauge fields, and hence can be pulled out of the path-integral Z K and 
absorbed in the overall normalization. The final result is then: 

Z K = N I VA exp [i ( Sym [A] + S^ [n ■ A] + JoA + K K [A]) } , (C14) 

where, from (C9), 



S&" [n ■ >t] = exp | -— / d"x [n ■ A a (x)]' f . (C15) 

In conclusion, the property of gauge field independence of the Fadeev-Popov determinant proves that there are indeed 
no ghost fields coupling to the gluon fields, hence the formulation is ghost-free. 
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APPENDIX D: GLUON PROPAGATOR AND POLARIZATION TENSOR IN THE AXIAL GAUGES, 

AND IN THE LIGHT-CONE GAUGE 



1. The general case 

In order to find the explicit form of the gluon propagator A K = Q K , we first evaluate its inverse, T K = (Q^ ) _1 , 
from the second functional derivative of T K with respect to A, and then invert it. With the conventions of Appendix 
B, we have, for the inverse of the exact gluon propagator A K , 



(A«)i (*>!/) = ^i 2 U 



S 2 T K [A] 



5 A (y)SA^(x) 



A=0 



2iST K [A = 0] 
SA VIJ .(y,x) 



(Dl) 



where the explicit form of the effective action r K [A] is given by (p4) and (p5J) together with the expressions rt5L)-(p3h. 
The exact propagator A K is related to the bare propagator A^ and the proper self-energy tensor H K through the 



Dyson- Schwinger equation 



(a^L = (^ - 1 ) + K„ 



n (o) 



n 



K, flV 



where, in the class of axial gauges, the propagator is transverse to the gauge vector n p 



TVAr = = AJfn,, 



(D2) 



(D3) 



while the polarization tensor is strictly transverse with respect to the external momentum q^ (the conjugate of <9 M ), 

3 M nr = o = ft, fir, (D4) 

and both are symmetric under interchange of indices and arguments, 

A£" = A v / , n^ = n^ . (D5) 



In order to infer the general form of the exact propagator A K , we apply (B7I) to (51 - 53), and carry out the 



Fourier transformation to momentum space. Then one observes that the axial-gauge representation of the inverse 
gluon propagator (D2), A^ 1 = H K + II K , can be decomposed into two independent Lorentz tensor components [p2[ , 



(A- 1 )^) = 5 ab (a K (q 2 , X )P^(q) + b K {q 2 , X ) Q^q)) 
with the projectors P^ v = P„ M and Q M „ = Q ufl , 



p ^{q) = 9W + 



w{q) = - 



1-X 



X— 2 — + X-^2~ 

n A <rt . rt r>^ 



n • q 



i-x 



q^q v _ n^ + q il n v 
q 2 n ■ q 



(1 - X )n 2 \ Uu,n 



i,'! 2 



VfA'VL/ 



(D6) 

(D7) 
(D8) 



which are orthogonal (in the space transverse to q) and obey the relations P^\P^ = P^v, Q^\T^ = Q^ v , P^xQi = 0. 
The invariant functions, a K and b K depend in general on q 2 and on the variable 



X 



X(n,q) 



n 2 q 2 
(n-q) 2 ' 



(D9) 



because the inverse propagator (D7) requires a scaling invariance under the change n^ -> An M . In a similar way, one 
may parametrize the propagator itself as 



VjUff) - 5 ah (A K (q 2 , X ) S^(q) + B K (q 2 , X ) T^{q)) 



(D10) 
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with different projection operators «S^„ = S ufJ _ and T^ v = T vli , 



S^(q) = 9u,u + — 
1 



T ^{q) 



i-x 



x (i + tf) Mil _ n » q » + ?/*"" + x "/*"" 

q 2 n ■ q n 2 



X l+§« j , 



n ■ q 



(Dll) 
(D12) 



which are again orthogonal (but now in the spa c e tra nsverse to g), satisfying S^xS^ = S M „, T^\T* — T llUl S^\T X = 0, 
and moreover n^S^ = = n^T^. Using (p7|)-(pl0|) in 



A Kj/lA (g) A K x „{q) = g^ , 



(D13) 



a^(q 2 ,x) ' 



S«(g 2 ,x) = X 



it is straightforward to derive 

A K (q 2 , X ) ■■ 
The bare propagator A„ corresponds to II = in (D2) and (D6), which yields a K = b K = q 2 + R K (q) 2 , and so, 

S ab 



M<7 2 ,x) 



(D14) 



a<°>° 6 (<z) 



q z + R K \q 2 ) 



§ ab 



while the inverse 



g 2 + iW) 



IIk reads 



n ii q v + q l xn v 2 2 g p g„ 



n • g 



(n-g) 2 



(D15) 



(a2») _1 

(A(°)^) Qb (g) = ^(g 2 + i? K (<Z 2 )) [ P M „ + 



6 ab (q 2 + R K (q 2 )) 



9/j.u 



q^qv 

q 2 



2 9 , J- \ TtuT^v 

n q ' ! 



e; (^-g) 2 



(D16) 



In order to derive expressions for the exact propagator A K and its inverse A K \ in correspondence to (D15) and 
(Dlq), it is useful inspect in more detail the structure of the polarization tensor II K , and which of its contributions 
are dominant. Let us define the dimensionless scalar functions 



nl 1} (g) = i - 



g K (q 2 ,x) 

q 2 + R K (q 2 ) 



4 2) (g) 



QK(g 2 ,x)-^(g 2 ,x) 
g 2 + i? K (g 2 ) 



(D17) 



and rewrite eq. (D6), on account of (D2) as 

iW<z,-g) = (a- 1 - Af)- 1 )^ = (a- 1 - 4°g 

= ( q, v q -^) n« + ( Ml _ M* + *»"* + (nV + r i } p^\ n ( 2 ) 

V 9 / \ f n-q (n-q) 2 J 

which implies the relations 

nw = i^ + ^V + r^T^^^nr n< 2 ) = 3nW-n^ 



{n-q)' 



(D18) 



(D19) 



In view of (D18), one realize s tha t II K consists of a covariant piece oc Tr K ' plus a non-covariant piece ex II« . 
Furthermore, comparing with (Dlq), in is obvious that it is solely the covariant contribution that survives in the limit 

of vanishing coupling g = 0, because then a K = b K = q 2 + R K , so that H K = 0. 



From ( D16 ) and ( D18 ), we read off the inverse of the exact gluon propagator (going over to £ — * 0), 
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(A- 1 )^) = 6"» ( q > + R K (<?)) { (gp, 



q»q» 



i-n 



2„2 



Qtj.qu n^q v + q^n v n*q* n li n v \ (2) 



and the actual gluon propagator is readily obtained by inverting (D20), 



K ab M 



8' 



(lb 



q 2 +R K (q 2 ) Vl-nl ] 



T (i) 



9V 



n ■ q 



n ■ q 



(n • qf 



2„2 



n^q v + q^riv q^q v I n z q 



n 



(2) 



q 2 \(n-q) 2 1 _ (II^ - II^ ) , 



n 



(D20) 



(D21) 



2. The case \ -» 



t(1) 



Inspection of the expressions ( D2C ) and ( D21 ) exhibit the relative importance of the contributions oc Tl K and 
Uk, ■ If the terms involving Ir K ' could be droppped, then both (D20) and (D21) would become simply the bare 



— 



counterparts (D15) and (Did) for £ — > 0, modulo the factors 1 — TL'k 1 , respectively 1/(1 — II« '). Now, there is no 

immediate argu ment why II^ itself should be negligable as com 
the propagator ( D21 ) could be small or even vanishing, is when 



T2i • • (1) (2) 

immediate argu ment why 11^ itself should be negligable as compared to H K , so that the only way the 11^ -term in 



2 9 

n^q* 
{n-qf 







which implies 



(n • q) 7 







(D22) 



(D23) 



The first condition corresponds to very large momentum component along the direction of n, for example, if n is 
chosen along the z-axis, then q z — > oo would do the job. The second condition, on the other hand, corresponds to 
picking, out of the class of axial gauge s, sp ecifically the light-cone gauge. Under either of these conditions, one arrives 
the very simple forms for ( |D20| ) and ( ]D21 ): 



(A^O?) = ^ (q 2 + R K (q 2 )) (1-Hx) ( g, v - ^) 



\dh 



K a M 



R K (q 2 



1-Hi 



9^ 



n ■ q 



(D24) 



(D25) 



One sees that now the effect of gluon self-interactions are encoded multiplicatively, so that we can express both A K and 
A^T 1 as the bare counterparts A^ , respectively A K , modulo a scalar renormalization function Z K (q 2 ,x)\x->Q = 
Z K (q 2 ), 

~ 2 ' R K (q 2 " 



(A-T"(a)- l (A^-- 1 Y b (q) 5 ab ( q2 + RK{q2) ) (q -^] 

{^)^W Z R (q 2 ) \*« )J q) - \ Z K (q 2 ) ) [^ q 2 ) 



(D26) 



A ab (a) 



M 2 ) U^X\q) 

\ /LIU 



tab 



Z K (q 2 



9i*v 



n^qv + qp,n v 



v*> \q 2 + Rn,{q 2 ) J V ' " l1, 9 

where the renormalization function Z K (q 2 ) is related to the gluon self-energy Tv K by 

Z K {q 2 ) 



1 



i - n«(g) 



with initial condition at q 2 = A 2 — > oo, 



Z «tf)\*= 



q 2 =A 2 



(D27) 



(D28) 



(D29) 



so to ensure that the full propagator equals the bare one in the asymptotic freedom regime when T1 K , H K — > 0. 
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APPENDIX E: SPECTRAL REPRESENTATION OF THE GLUON PROPAGATOR IN THE AXIAL 

GAUGES 



In this Appendix we discuss in more detail the relation between the gluon propagator A K and its spectral density 
p K , a s intr oduced in Sec. 3, eq. (110). Recall, that the gluon propagator is formally defined, according to (B2), (Bll) 
and (B18), as the connected 2-point Green function in the presence of the infrared cut-off k, involving the time-ordered 
product of two gauge fields at space-time points x and y Fj: 



A »*(*,*/) = (0|T \A a Jx)Al(y) |0), 



(El) 



Analogously we define now the gluon correlation function as the non-time ordered 2-point function which describes 
the correlation between two gluon fields at x and y, irrespective of their time history and spatial origin: 



P*Z,(x,v) = <0| [A a Jx)Al(y)] \0) K 



(E2) 

(E3) 
(E4) 

We recall that both the propagator and the correlator depend on q and n, more precisely on q 2 and n ■ q. Let us 
now focus on the correlat ion function /j Ki/jj , and then work our way back to the propagator A K-/J „. Following pl|, we 
take the commutator in (E2) apart and define (suppressing the color indices, as they are in paralell with the Lorentz 
indices) , 



In momentum space, we write (El) and (E2) as, 

A K ,nu(q) = fd 4 xe tqx A K ,^(x,0) 
Pk,^u{q) = d 4 xe iqx p K ,p U (x,0). 



p { kXM) 



Ae"' 1 (0|^(x)X(0)|0)« 



d 4 xe lq - x (0| A(0)^(x)|0) 



Hence, 



PK,ixu{q) = P^Lil) + pi.L(l) 



dW-(0| [A^(x)A b u (0)] |0) K , 



and we have the following properties, 



PK,nv{q) = -pK tlJ ,u{-q) 



(E5) 
(E6) 

(E7) 
(E8) 



Pk,Iu(q) = -p^lui-i) 



(E9) 



Now, recall that in the axial gauges the g-dependence of both the propagator A K and the correlator p K can enter 
only in terms of the two invariants q 2 and (n • q) 2 . It is therefore useful to introduce a notation for the decomposition 
of an arbitrary four- vector v^ into its longitudinal (v^) and transverse components (i>j.) with respect to the gauge 
vector 7J M : 



(n ■ v) n^ , 



(E10) 



with v\ = (n ■ v) n 2 , v 2 
this notation, 



v\, and n ■ vl = n ■ vt = 0. Thus, the g-argument in p K , for instance, reads with 



L We suppress here the superscript (c) for 'connected' 
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PK,tiu{q) = pK,nv(q ,q L ) ■ 



(Ell) 



In order to derive the relati on b etween the time-ordered product of gauge fields (El) in the propagator A K and 
the non-time-ordered product (E2) in the correlation function p K , we proceed now as follows. Let {| N)}^ denote 
a complete set of states which spans the Hilbert space Tie of all possible gluon configurations (one of which is the 
vacuum state | ) ) . Inserting then this complete set of gluon states into ( E5 ) gives 



P^iM = E <0M M |iV> re <iV|A,|0> K {2n) i 8\q-p N ) 



(E12) 



and require qo — po > 0. Inserting ( E12 ) into (E5) and inverting the Fourier transform, one readily finds, 

d 4 q 



;o|^(x)A(o)|o) K 



(2tt)4 

A 

(2tt)* 



e- 1 *" %o) p K ,^{q\q 2 L ) 



* %o) / dy L e'™ p K ,^(q 2 ,y L ) 



(E13) 



Ak , respectively, 



where p K (g 2 ,j/L) is the longitudinal transform of ( |E11 ). Introducing the advanced and retarded functions, A K and 

Ai±)(x) = ± J -0-c c-*«- %o) 5(q 2 - n 2 ) , (E14) 



one can express ( E13 ) as 

(0|^ M (a;)^(0)|0)« 



dq 2 / dy L p K ,^(q 2 ,y L ) A^ +) (a; - y L ) 



Similarly, from the crossing relations (E8) and (E9), one obtains for the reversed product of the gauge fields, 



;0| A v {fi)A li {x)\Q) K 



dq 2 



dyL p K ,^u{q 2 ,yL) A£ \-x + y L ) 



(E15) 



(E16) 



With the above relations we can now express the time-ordered product of the gauge fields, which determines the 
propagator A K via (El), as PI 



;0|T [AXx)Al(0)] |0) K = 



dq 2 



dy L p K ,n V {q 2 ,y L ) \ 9{x Q ) A[+\x - y L ) + 9{-x ) A[ \x-y L ) 



(E17) 



If the gauge vector is chos en to be space-like or light-like, i.e., n 2 < 0, then causality allows us to replace 9(xq) by 
0(xq — y^, in which case (E17) together with (El) yields 



A K>ftv (x,0) 



dq 2 / dy L p K ,^(q 2 ,yL) A^(x-y L ) 



An 



where A K is the standard Feynman function, 

A[ F \x) = 9(x )Ai+\x)+ 9(-x Q ) A[r\x) . 



(E18) 



(E19) 



12 To be precise, here the indices i,j = 1,2,3 should be restricted to the sp atial components of the gauge fields A and of 
p re , because in the coordinate representation, the tensor structure of p K , M ^ (cf. (E22J below) leads to space-time derivatives d^ 
acting on the A„ ' functions, which causes the time-ordering operation does not commute with the time derivatives arising 
from the time components of p Ktlxv . 
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In momentum space, eq. (E18) reduces to the well-known spectral (or Lehmann) representation, 

„/2 Pn,^(q 2 ,q 2 L ) 



A Kl/iI/ (g) = / 
Jo 



dq' 



1 i ? 

q z — q ' z 



If we decompose the tensor structure of A K ^ as in Sec. 3, eq. (|85|), or Appendix D, eq. (D10), 

A„ l/U ,(«) = A K (g 2 , X ) V(?) + B K {q\ X )T^{q), 
with x = n 2 q 2 /(n ■ q) 2 as before, and analogously, for the correlation function p Kt/ j,v, 

PK,v,»{q) = pi(q 2 ,X) S^{q) + p^(q 2 ,x)T flu (q) , 



then eq. ( E20| ) may be written as (noting that \ = n 2 q 2 /{n ■ q) 2 = n 2 q 2 jq 2 Ll i.e. \ oc g L 2 ), 

,2Pn(q 2 ,x) 



A-u^(g) = 5m, (g) 



q A — q l 



dq 



q 2 — q' 2 



(E20) 



(E21) 



(E22) 



(E23) 



In the case x - ► 0, corresponding to q\ — > oo, the second term in (E23) tends to zero, as discussed in Appendix D. 
Thus, with 



pi{q\x) 



P.(q 2 ) 



P B K {q\x) x ^ o 



we are left with 



A K ,^(g) = S'{q) / dq 



/2 P«(9' 2 ) 



(E24) 



(E25) 



which is precisely eq. (|110|) with the substitution q 



R K in the denominator. 



The spectral representation ( E23| ) or (E25) has a rather intuitive physics interpretation: The propagator for a gluon 
of momentum q is a sum over all intermediate virtual gluon states of momentum q' with virtuality q 1 2 and with level 
density p K {q'), and weighted with the phase-space factor l/(<7 2 — q' 2 . For a bare gluon with momentum q — > oo in 
the ultraviolet, no virtual fluctuations would be present, so that p K (q' 2 ) = S(q' 2 ). On the other extreme, an infrared 
gluon with momentum q — > is dressed by a dense cloud of soft virtual fluctuations, so that p K (ci' 2 ) can be a very 
broad distribution. 
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